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ABSTRACT 


A new concept in nuclear power generation is being explored 
which essentially eliminates major objections to nuclear power. 

The Satellite Nuclear Power Station (SNPS), remotely operated in 
synchronous orbit, would transmit power safely to the ground-by a 
microwave beam. Fuel reprocessing would take place in space and 
no radioactive materials would ever be returned to earth. Even 
the worst possible accident to such a plant should have negligible 
effect on the earth. 

An exploratory study of an SNPS power station to provide 
10,000 MWe to the earth has shown that the system would weigh 
about 20 million pounds and cost less than $1000/KWe. An advanced 
breeder reactor operating with an MHD power cycle could achieve an 
efficiency of about 50% with a 1100°K radiator temperature. If a 
hydrogen moderated gas core reactor is used, its breeding ratio of 
1.10 would result in a fuel doubling time of a few years. A colloid- 
core or NERVA type reactor could also be used. The efficiency of 
power transmission from synchronous orbit would range from 70% to 80%. 
The only environmental effect of this power plant would be a slight 
thermal discharge at the receiving antenna, equal to about 10% of 
the heat released by today's most efficient power plants. Thus, the 
SNPS comes close to the ideal of economical power without pollution. 



/Yv 3 - 1 . y 


3 


INTRODUCTION 

J * R. Williams 


Preceding page blank 



5 




Preceding page blank 

WHY SNPS? 

It is generally acknowledged today that the world is 
facing an energy crisis. Electric power requirements have been 
doubling every ten years, and the demand for transportation and 
heating fuels has been increasing rapidly. The rapidly increas- 
ing demand for fossil fuels has pushed prices up and reduced 
their availability. In view of the higher costs of fossil fuels 
and increasingly tight restrictions on the emission of atmos- 
pheric pollutants from fossil fuel combustion, nuclear power 
has become competative and hundreds of nuclear plants are built, 
planned or under construction. However, grave questions are now 
being raised about the safety and possible adverse environmental 
effects of these power plants. 

The main objections to nuclear power are as follows: 

(1) Radioactive emissions during normal plant operation - This 
is real ly not a problem in that these emissions can be reduced 
so that exposure to the public is far below background. This 
objection can be resolved by proper plant construction, proper 
plant maintenance, and public education. (2) Thermal pollution 
from nuclear power plants - This is a problem with all thermal 
power plants (including geothermal). Heat rejection to rivers 
upsets the ecology of the rivers, wet cooling towers produce 
local fogging conditions, and dry cooling towers produce thermal 
plumes which can be a significant hazard to aircraft and which 
also effect local meteorological conditions. The ultimate heat 
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rejection method may be radiation to space. However, as the 
heat rejection temperature is raised, plant efficiency drops 
and the total thermal discharge increases. Various ideas have 
been advanced about how the waste heat might be used constructively, 
but it is difficult to find a practical use for such large amounts 
of low grade heat. (3) Accidents involving a reactor or a fuel 

reprocessing plant which result in a release of fission products 

are a major concern. If such an accident occurs, it could deal 
a severe setback to the development of nuclear power. (4) Fuel 
element shipping is a major problem with respect to nuclear power 
plants. Highly radioactive fuel elements must be removed from the 
power reactors and transported to a reprocessing plant. An accident 
which released fission products or plutonium would constitute a 
major hazard to the public. (5) Safeguards present perhaps 
the greatest long-term problem of nuclear power plants. The 
amount of plutonium in a single fast breeder reactor is suffi- 
cient for the construction of over a hundred atomic bombs. If 
these power plants are to be a major worldwide energy source 
for the future, how can one insure that a small nation or group 

of individuals does not divert any plutonium to the production 

of weapons? This is a problem that may be extremely difficult 
to resolve. The general availability of plutonium could be a 
major factor in the proliferation of nuclear weapons. 

(6) Disposal of radioactive wastes is presently of major concern 
to environmentalists. How can one be sure that buried wastes 
will, at no time in the future, leak radioactivity into the 
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environment? In spite of the elaborate safety precautions that 
are taken, local governments tend to be strongly opposed to the 
disposal of radioactive wastes in their area. One major reason 
alternative energy sources are being persued with such vigor is 
to reduce or eliminate radioactive waste disposal on the earth. 

A new concept in nuclear power generation is being explored 
which essentially eliminates all six objections to nuclear power. 
The Satellite Nuclear Power Station (SNPS) would be located in 
synchronous orbit and transmit power safely to the ground by a 
microwave beam. Fuel would be reprocessed at the plant or 
elsewhere in space, and no radioactive materials would ever be 
returned to the ground. Even the worst possible accident to 
such a plant would have negligible effect on the earth. Radio- 
active wastes would be placed in containers and dumped into 
the sun or placed in a solar orbit inside the earth's orbit. The 
safeguards problem, which is a very difficult problem for future 
ground based nuclear plants , is greatly reduced since the SNPS 
plant would be in synchronous orbit and could be remotely operated. 
The resolution of this and the other drawbacks of nuclear power 
may be well worth the additional cost of power from the SNPS. 

If the SNPS is to serve as a major energy source for the 
future, it must be a breeder so as not to deplete available fuel 
resources. The plant should have a reasonable thermal efficiency 
even though heat must be rejected by a radiator. It must be 
large, of the order of 10,000 MWe or more, so that the unit cost 
of delivered electric power is reasonable. The fuel cycle must 
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be as simple as possible to reduce the cost and complexity of 
on-site fuel reprocessing. A nuclear power plant which has 
been under study and which appears to meet these requirements 
would use an advanced high-temperature breeder reactor with MHD 
energy conversion. 

In view of the current uncertainty in the future avail- 
ability of fusion power and the difficulties facing terrestrial 
fission power plants, the SNPS should be considered a major 
alternative energy system for supplying tommorrow's energy needs. 
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ADVANCED REACTORS FOR NUCLEAR-MHD 

In order for efficient multi -megawatt closed cycle 
nuclear-MHD systems to become practical, long-life gas cooled 
reactors with exit temperatures of about 250Q°K or higher 
must be developed. Four types of nuclear reactors which have 
the potential of achieving this goal are the NERVA-type solid 
core reactor, the colloid core (rotating fluidized bed) reactor, 
the "light bulb" gas core reactor, and the "coaxial flow" gas 
core reactor. 

1 ? 

The solid core NERVA type reactor, * which is already well 
developed, offers the promise of almost immediate application for 
MHD power generation. The colloid core reactor^ has been 
studied by the Air Force Aerospace Research Laboratories for the 
past eight years, and their developmental program has now reached 
the point that a contract has been given to the Battel le Memo- 
rial Institute for an in- reactor test of a fission-heated colloid 
core reactor experiment using UOg particles in a confined vortex. ^ 
This two year experimental study is the logical step prior to 
the development of a full scale colloid core reactor. The 
colloid core reactor uses a rotating fluidized bed of uranium 
dioxide particles in a confined vortex to heat a gaseous working 
fluid to as high as 3200°K, temperatures which are ideal for 
closed cycle MHD power generation. The nuclear fuel cycle, in 
comparison with present fuel cycles, is greatly simplified since 
there are no fuel elements. 
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Two types of gaseous core nuclear reactors also show 
promise for MHD power generation, the nuclear "light bulb" 
reactor 6-10 and the coaxial flow reactor.** -11 In the light 
bulb reactor, gaseous nuclear fuel is confined within a trans- 
parent partition and the working fluid is heated by the absorp- 
tion of thermal radiation transmitted through the transparent 
partition from the fissioning gaseous fuel. 

Prior to the recent NASA cutback in January, 1973, the 
United Aircraft Research Laboratories was preparing to conduct 
a small scale fission heated light bulb reactor experiment in the 
Nuclear Furnace reactor. 12 Uranium gas was to be confined in 
a transparent partition and heated by fission to a very high 
temperature, while a gas such as argon, helium, or hydrogen 
flowing around the partition is heated to about 3500°K by the 
thermal radiation from the hot uranium gas inside the partition. 
The NASA-Lewis Research Center was also proceeding with plans for 
a Fissioning Uranium Plasma Test Facility to be located at the 
Nuclear Rocket Development Station. This reactor would have 
used MTR type fuel elements surrounding a two-foot diameter 
cavity, and was to be used to test the various gas core and 
colloid core systems, to demonstrate MHD power generation with 
these reactors, and to study other applications of fissioning 
uranium plasmas. Since it now appears that NASA will no longer 
be involved in the development of nuclear reactors, it is hoped 
that another agency will continue the development of these 
high temperature reactor systems for power generation. If such 
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development is continued, construction of prototype test reactors 
could begin very soon and they would probably be operating by 
1980. The fissioning uranium plasma test facility proposed by 
NASA would cost about 16 million dollars to build, and could be 
used to confirm the technical feasibility of larger colloid core 
and gaseous reactor systems, and to study the performance of MHD 
generators operating with these reactors. 



SNPS POWER PLANT SYSTEMS 


For the past 1 1/2 years the authors, with NASA support, 
have been evaluating MHD power plant systems utilizing these 
high temperature reactors. Some preliminary results have been 
reported, 13,14 anc j earlier papers^-2l describe previous 
studies of gas core reactor MHD power plant concepts. Some of 
these earlier studies^*^ were aimed at determining whether 
or not a gas core reactor can breed its own fuel. The first 
calculations considered gas core fast breeder reactors, and 
showed that although the breeding ratio was high, the critical 
mass was also large. Gas core thermal breeder reactors, 
moderated by hydrogen gas, were shown to have much lower critical 
masses and reasonable breeding ratios. 

Three different types of closed cycle nuclear MHD power 
plant systems have been analyzed to determine the operating 
characteristics, critical parameters, and performance of these 
power plant systems. The basic power cycles which have been 
studied are illustrated by Figures 1-3. Each of these power plant 
systems may be subdivided into three component subsystems 
(Figure 4): 1) the high temperature reactor with attached MHD 

generator and uranium separator (if required), 2) the 
compressor system and 3) the heat rejection system, which 
is a radiator. 

The first subsystem, which is the same for all 3 plant 
configurations studied, consists of the nuclear reactor, the 
MHD generator, uranium separator (if required) and all associated 
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MODE I 


Figure 1 


Turbine-Compressor Cycle with High Temperature Regenerator 
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MODE II 

Figure 2. Turbine-Compressor Cycle Without High Temperature Regenerator. 













MODE III 

Figure 3. Motor-Compressor Cycle With High Temperature Regenerator 
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Electric power produced by 
the MHD generator. Subsystem 
contains all nuclear components 


Electric power generated or . 
required, depending on cycle 
and operating conditions. 


Some electric power required 
by liquid metal pumps. 



Figure 4. The Three Major Power Plant Subsystems Common to All 
Three Cycles 



Figure 5. Four Pass Tube-Fin Gas-to-Sodium Crossflow Heat 
Exchanger 
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uranium recycling and reprocessing facilities. Figures 1-3 
show a gas core reactor of the "coaxial flow" type, but any of the 
four reactor systems could be used. Both the coaxial flow gas 
core reactor and the colloid core reactor would require uranium 
separators, as shown in figures 1-3. The light bulb and NERVA 
type reactors would not require separators, since the uranium 
fuel would not become mixed with the working fluid. This first 
subsystem contains all the nuclear components of the power plant 
and the components that require the most technological develop- 
ment. These are the components that would be developed with a 
uranium plasma test facility of the type proposed by NASA. 

The second subsystem consists of the turbine, compressor, 
and associated heat exchangers for Modes I and II; and the 
compressor, electric motor and heat exchangers in the case of 
Mode III. In Mode I (figure 1) a regenerative heat exchanger 
is used to cool the gas from the MHD exit temperature to an 
acceptable turbine inlet temperature, while the compressed 
gas returning to the reactor is heated. In Mode II, (figure 2) 
the regenerator is removed and the temperature of the gas 
exiting the MHD generator is reduced to the turbine inlet temp- 
erature by mixing with cooler gas from the first stage compressor. 
This avoids the problems associated with the high temperature 
regenerator, but at the expense of cycle efficiency. Cooling 
is provided by gas-to-sodium tube-fin heat exchangers. Mode III 
(figure 3) uses a high temperature regenerator but eliminates the 
turbine. The major advantage of this cycle is that there are no 
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moving parts at high temperature, and the efficiency is only 
two or three percent less than Mode I. 

In general, the Mode III cycle appears to be the most 
attractive because of its simplicity and potential for high 
reliability, but it will require the development of efficient 
high power (probably cryogenic) electric motors. Mode I is the 
most attractive cycle if such motors are not developed, and 
provides the highest cycle efficiency. However, if regenerator 
problems prove insurmountable. Mode II can be used. Mode II 
can use current technology components for this subsystem. 

The third subsystem rejects the heat removed by the liquid 
sodium from the sodium heat exchangers. In space the sodium is 
circulated through a heat-pipe radiator and the heat is rejected 
to space. 

Figures 1-3 show a "coaxial flow" gas core reactor in the 
first subsystem, although any of the other three reactor types 
could be used. Uranium fuel separators would not be needed with 
the solid core or the "light bulb" gas core reactor. All these 
reactors, except the solid core, require continuous fuel recircu- 
lation, and also permit continuous fuel reprocessing and the 
removal of gaseous fission products. The probable reactor 
operating temperature range is given in Table 1. 
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Table 1. Reactor Exit Temperatures for Advanced Power Reactors 


Reactor 


Temperature °K 


Solid Core (NERVA type) 
Colloid Core 
"Light Bulb" Gas Core 
"Coaxial Flow" Gas Core 


2200°K - 2500°K 
3000°K - 3200°K 
3500°K - 4000°K 
3700°K - 5000°K 


The compressor subsystem uses either a turbine (Modes I and 
II) or a cryogenic electric motor (Mode III) to drive the multi- 
stage compressor. Cylindrical plate-fin counterflow surface 
compact heat exchangers are used for regeneration and four pass 
gas-to-sodium crossflow type heat exchangers (Figure 5) are used 
for primary heat rejection and intercooling between compressor 
stages. The surface characteristics of these heat exchangers 
are given by Tables 2 and 3. 


Table 2. Gas-Gas Heat Exchanger Characteristics 


Hot Side Cold Side 


Surface 

plate- 

fins 

plate- 

fins 

Plate spacing (ft) 

0.25 

0.204 

Hydraulic radius (ft) 

0.00253 

0.000943 

Fin thickness (in) 

0.006 

0.006 

Heat transfer area/vol . 
(ft 2 /ft 3 ) 

367 

855.6 

Fin area/total area 

0.756 

0.884 






Figure 6. Mode III SNPS Power Plant 













This page is reproduced at the 
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reproduction method' to provide 
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Figure 7. Satellite Nuclear Power Station in Synchronous Orbit. 
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Table 3. Gas-Na Heat Exchanger Characteristics 



Gas Side 

Na Side 

Surface 

conti nuous 

fl at 


fin 

tubes 

2 

Frontal Area per tube (in ) 


0.434 

Fin thickness (in) 

0.004 


Free flow area/frontal area 

0.780 

0.129 

Fin area/total area 

0.845 


Hydraulic radius (ft) 

0.00288 

0.00306 

Heat transfer area/vol. 

(ft 2 /ft3) 

270 

42 


Experimentally determined correlations between Reynolds number 
and friction factor and heat transfer characteristics are used 
to evaluate the performance of the heat exchangers for each 
specific plant operating condition. 

An artist's concept of a MODE III SNPS power plant is 
illustrated by figure 6. Figure 7 depicts this power plant 
in synchronous orbit with microwave transmission to a receiving 
antenna on the earth. Studies at Reytheon^ and the Grumman 
Aerospace Corporation"^ have shown that safe microwave transmission 
of electric power from synchronous orbit is feasible with efficien- 
cies of 70 to 80%. These studies have been performed in connection 
with a Satellite Solar Power Station (SSPS) study currently under- 
way. The SNPS would have the advantage of SSPS in providing 
power without pollution. 
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MHD GENERATOR PERFORMANCE LIMITATIONS 


R. J. Rosa 
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MHD GENERATOR PERFORMANCE LIMITATIONS 

The electrical properties of the gas are the primary 
determinant of whether a generator can be operated at a high 
loading factor without excessive length or field strength 
requirements. The electrical properties of a plasma that are 
relevant to MHD are the conductivity and the Hall parameter. 
Assuming that a given level of electrical power is sought, and 
that the generator L/D is fixed, for each pressure there is a 
minimum temperature that is necessary. The parameter L/D is 
determined largely by the boundary layer growth. Experience 
has shown L/D-10 to be about right in inert gas generators. 

Figures 8 and 9 assume that L/D=10, and show the maximum 
allowable pressure at several given levels of power extraction, 
vs. temperature. These plots show the basic advantage of higher 
temperature as it relates to conductivity; there are two 
additional factors to be considered. 

First, as temperature increases we can increase pressure to a 
level higher than before available, and boundary layers become 
much better behaved. L/D = 10 is probably a conservative estimate 
of what is allowable. 

Secondly, the higher pressures available above -3500 °K lower 
the Hall parameter, «r, so that a continuous electrode generator 
becomes a possibility. The power extracted from a conti nous 
channel as opposed to an infinitely finely segmented one is given 
by the factor 1 

1 + (tUTJ 



(atm) 
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T 01 (thousands of °K) 


Figure 8. Maximum Allowable Reactor Pressure and Minimum 
Temperature to Preserve an Acceptable Duct L/D 
for a Pressure Ratio of 30. 





Max. Allowable Reactor Pressure 
and Min. Temperature to Preserve 
an Acceptable Duct L/D 

H 2 ♦ 1% Cs or Equivalent 
Assumed Molecular Wt. * 4 
B = 10 

3 


1 


1 
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where the Hall parameter, uyj-t is proportional to the magnetic 
field strength and inversely proportional to pressure. For a 
reactor temperature of 3500°K and pressure of 200 atmospheres , 
ccris about 0.1 at the inlet and 0.8 at the exit. This implies 
an average power differance of only about 10%, so the continuous 
electrode generator would be feasible. For higher reactor 
temperatures the difference would be even smaller. 

The use of continuous electrodes would eliminate the worry 
about the electrical integrity of the electrode design. This 
would be of considerable practical importance. 
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NUCLEAR ANALYSIS 


J. D. Clement and K. D. Kirby 
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SUMMARY OF NUCLEAR CALCULATIONS 

Exploratory calculations have been performed for several 
gas core breeder reactor configurations. The computational 
method involved the use of the MACH-1 one-dimensional diffusion 
theory code and the THERMOS integral transport theory code 
for thermal cross sections. Computations have been performed 

to analyze thermal breeder concepts and non-breeder concepts. 

?33 

Analysis of breeders has been restricted to the U -Th 
breeding cycle, and computations have been performed to examine 
a range of parameters. These parameters include U^ to 
hydrogen atom ratio in the gaseous cavity, carbon to thorium 
atom ratio in the breeding blanket, cavity size, and blanket 
size. Results of a parametric survey show that breeding ratios 
in the range of 1.06-1.12 could be obtained with critical masses 
of 300 to 850 kilograms U233 for various material compositions in 
a 5 meter diameter cavity with a 0.5 meter thick blanket. The 
effect of fissile material in the blanket, cavity temperature, 
and structural material in the blanket has been estimated. The 
breeding ratio can be increased to 1.13 by utilizing fissionable 
material in the blanket without a large increase in total U 233 
mass. A decrease in average cavity temperature from 4000°K 
to 3000°K increases the breeding ratio from 1.10 to 1.12 with a 
significant reduction in cavity pressure. Cavity pressure at 
3000°K is about 400 atmospheres. Structural material decreases 
the breeding ratio by approximately 2 % for 0.2 atom percent 
natural molybdenum or 4% enriched molybdenum in the blanket. 
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Gaseous core reactors, non-breeding in nature, were also 
analyzed with different fuels and for varying sizes. Cavity 
diameters ranging from 1.2 to 3.0 meters with BeO reflectors 
0.3 and 0.5 meters thick were examined with U 233 fuel and U 235 
fuel of various enrichments. Results show U critical masses 

por 

significantly lower than U critical masses due to the low 
energy fission resonances in However, for high enrichment 

. OOC 

(*9335) the U requirements are less than 15 kilograms. 

Pressure for the larger cavity sizes is generally below 300 
atmospheres for U 233 or highly enriched U^. 
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COMPUTATIONAL METHODS 

Nuclear analysis of the gaseous core nuclear reactor is 
a very difficult task requiring highly sophisticated techniques. 
Several analyses C *~ CK> have been performed which have used 
very sophisticated techniques and pointed out the areas of 
difficulties. For examining a broad range of designs however 
one may utilize less sophisticated techniques to observe trends 
and perform parametric studies in order to identify concepts for 
further study. 

The first phase in performing exploratory nuclear analysis 

for the gaseous core nuclear reactor involved implementing the 

necessary computational tools and formalizing a computational 

method. The major portion of the effort early in the study was 

devoted to this area. In order to expedite this phase the MACH-1 
27 

code was used as the primary computational tool in the nuclear 

analysis. To allow a more realistic model of thermal neutron 

processes in the high temperature gaseous core reactor concept, 

28 

the THERMOS code was implemented to supply thermal neutron 
parameters to MACH-1. 

The computation method used in the nuclear analysis relies 
on these two codes. MACH-1 is a one-dimensional diffusion 
theory code with one thermal group (no upscatter) and THERMOS 
is a one-dimensional integral transport theory code in the 
thermal range with complete upscatteri ng. All reactor con- 
figurations are assumed to be spherical and hence amenable to 
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one-dimensional analysis. For the MACH-1 code the 26-group 
"ABBN" cross section set of Bondarenko, et al? 9 was used. The 
thermal group of the ABBN set is for 2200 m/s (0.0253 eV) neutrons 
and hence is not realistic for the high temperatures of a gaseous 
core reactor (5000°K, kT=0.43 eV). The THERMOS code was thus 
used to determine thermal cross sections to be inserted into the 
MACH-1 computation along with the ABBN set. For a given con- 
figuration the computational method was as follows: 

1. Run MACH-1 with 26-group ABBN to estimate critical 
concentrations and preliminary results. 

2. Run THERMOS with 50 groups (up to 2.15 eV) using 
above concentrations and calculate spatial and spectrum 
averaged cross sections. 

3. Run MACH-1 with 22 fast groups from ABBN (>2.15 eV) 
and thermal cross sections from THERMOS run. 

Thus the final results of a computation may be thought of as 
a 23-group calculation with one thermal group using a thermal 
cutoff of 2.15 eV. A schematic of the computational method 
is shown in Figure 10. Steps 2 and 3 could be repeated if final 
concentrations vary markedly from the estimates; steps 1 and 2 
could possibly be omitted for very similar configurations. The 
high thermal cutoff value is required because of the possibility 
of a large increase in neutron energy due to upscatter from the 
high temperature hydrogen moderator/coolant. 

Explicit in all calculations are the assumptions associated 
with the two computer codes. Diffusion theory does not seem to 
be very restrictive based on previous comparisons to transport 
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Step 1. Run MACH-1 with 26 group 
ABBN cross section set to esti- 
mate critical concentrations and 
get preliminary results. 


Step 2 . Run THERMOS with 50 groups 
(up to 2.15 eV) using above con- 
centrations and calculate spatial 
and spectrum averaged cross sections. 


Step 3 . Run MACH-1 with 22 fast 
groups from ABBN (>2.15 eV) and 
thermal cross sections from the 
THERMOS calculations. 


Figure 10 


Computational Method for Nuclear Analysis 
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theory for a fast reactor configuration (k correct! on=+. 009) . 
THERMOS contains the assumption of isotropic scattering but 
this is felt to be quite sufficient at the energies involved 
(<2.15 eV). More restrictive assumptions for the THERMOS runs 
are probably the si owing-down source and the 33 resonance 
below 2.15 eV. 

The si owing-down source for THERMOS is assumed to be 
spatially independent, MACH estimates as shown in Figure 11 show 
that the epithermal flux is rather flat in the cavity but 
decreases rapidly in the blanket region. This would imply 
then that the flat source assumption is rather good for the 
cavity and perhaps not as good in the blanket. But since the 
temperature is not as high in the blanket and resonance capture 
is important in thorium, results should not be as sensitive to 
thermal cross sections for blanket materials as for the cavity 
material . 

The THERMOS Code must also handle the resonances at 
1.78 and 1.55 eV since they lie below the thermal cutoff. No 
Doppler broadening capabilities exist with the code so these 
resonances are included at room temperature only. These indirect 
assumptions of no Doppler broadening of these resonances should 
not be too severe since the resonances are very broad even at 
room temperature. Since only eight of the fifty THERMOS groups 
are used to span these resonances, results are probably less 
sensitive to Doppler broadening than to the low number of groups 
in that interval. 
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Figure 11. Spatial Distribution of Epithermal Flux (E>2.15 eV) 
for a Typical Gaseous Core Breeder Confuguration 
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A general assessment of the computational method and its 

assumptions was provided by a check calculation on a configuration 

24 

analyzed by Whitmarsh. This case is the 10 ft. cavity diameter, 

2 ft. reflector region configuration described in Reference 24. 

An essentially equivalent configuration was obtained by reducing 
the number of regions by homogenizing similar regions. Then the 
computational method outlined previously with MACH-1 and THERMOS 
was used to analyze this configuration. The THERMOS computation 
was performed for the cavity regions only. Final results gave a 
value of k=0.986 for this configuration. In light of the homo- 
genization used to obtain a nearly equivalent configuration, the 
agreement tends to show the computational model to be valid. The 
largest source of discrepancy was attributed to the sensitivity to 
U235 thermal cross sections. This points out the need for a multi- 
thermal group treatment. The agreement does show that this 
computational method should be sufficient to identify trends and 
perform parametric studies for various gaseous core nuclear reactors. 
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GAS CORE BREEDER CALCULATIONS 

In this section results of the nuclear analysis of several 
concepts of a gaseous core breeder reactor are given. The 
primary objective of this portion of the nuclear analysis has 
been to perform nuclear calculations on various reactor con- 
figurations to determine a feasible gaseous core, thermal 
breeder, reactor power plant. Only thermal breeder configu- 
rations based on the Jh232_y233 breeding cycle have been 
examined. Although a fast breeder reactor may yield a higher 
breeding ratio, as found from a preliminary survey, the thermal 
breeder has the advantage of a much lower critical mass, simpler 
control, and in general, lower cavity pressure. If one uses 
the reactor doubling time (time necessary for the excess fuel 
bred to equal a new critical loading) as the figure of merit, 
the thermal breeder can compete favorably with the fast breeder. 
(The doubling time is directly preportional to the critical mass 
and inversely proportional to the breeding ratio minus one). For 
an extraterrestrial plant where excess fuel is desired a low 
doubling time is desired, but if all that is desired is to keep 
the original plant operating, then a larger doubling time (lower 
breeding ratio) merely compensating for process losses would be 
suffi cient. 

Since the thermal breeder does appear to be able to compete 
with the fast breeder and has advantages which could allow easier 
adaptation for extra-terrestrial use, nuclear analysis of several 
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configurations was carried out. 

As noted previously, all the configurations examined were 
spherical in geometry. These cases were described as two or 
three region spheres in the MACH runs and as slabs in the THERMOS 
calculations. The configurations examined are basically as that 
depicted in Figure 12. The cavity region contains hydrogen as 
moderator/coolant, U 233 as fuel, and sometimes thorium as fertile 
material. The blanket consists of graphite and thorium. The 
relative concentrations of the materials as well as the size of 
the regions were varied parametrically to examine a matrix of 
cases in an attempt to obtain the most feasible gaseous core, 
thermal breeder comcept. 

The first two parameters examined were the hydrogen to 
uranium atom ratio in the cavity (H/U ) and the carbon to thorium 
atom ratio in the blanket (C/Th). Initially cases were to be 
examined with H/U ratios ranging from 40/1 to 140/1 and C/Th 
ratios ranging from 2/1 to 50/1. Carbon atom density was kept 
constant in all calculations. Step 1 in the computational 
method (MACH-1 estimates) suggested that C/Th ratios greater than 
10/1 yield very low breeding ratios, and that H/U ratios below 
60/1 were undermoderated, hence the combined calculation (MACH-1/ 
THERMOS) was performed for H/U from 60/1 to 140/1 and C/Th from 
2/1 to 10/1. For all cases the cavity radius is 250 cm and blanket 
thickness is 50 cm. Results of these calculations are shown in 
Table 4 and Figures 13-15 for the important parameters of reactor 
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Figure 12. Typical Reactor Configuration for Nuclear Analysis 
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Table 4 Parametric Study of Relative Material 
Concentrations in a Gaseous Core Breeder Reactor 


Cavity 
H/U Ratio 

Blanket 
C/Th Ratio 

U 233 Mass 
(Kg) 

Breeding 

Ratio 

H Press 
(atm) 

Doubling ^ 
Time (yr) 

140/1 

2/1 

452 

1.1026 

(a) (b) 
710 514 

9.6 


4/1 

390 

1.0962 

612 443 

8.9 


10/1 

301 

1.0636 

472 342 

10.3 

100/1 

2/1 

576 

1.1056 

646 468 

11.9 


4/1 

494 

1.0997 

553 401 

10.8 


10/1 

375 

1.0662 

420 304 

12.4 

60/1 

2/1 

847 

1.1029 

569 413 

17.8 


4/1 

721 

1.0966 

485 351 

16.2 


10/1 

537 

1.0635 

361 261 

18.4 


(Cavity radius 

250 cm. Blanket 

thickness 

50 cm) 



(a) Hydrogen partial pressure at 4000°K, H„ mole fraction 0.92. 

(b) Hydrogen partial pressure at 3000°K, \\i mole fraction 0.99. 

(c) For 1000 Mw(t), proportionally lower per higher average power. 


breeding ratio 


3 


1.20 


C/Th = 2/1 


1.05 


C/Th 

'« 4/1 

C/Th= 10/1 

1 


LOO 


H/U atom ratio 


Figure 13. Effect of Thorium Concentration in the Blanket 
on Breeding Ratio 
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Pressure (atm) 
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Figure 15. Effect of Hydrogen/Uranium Atom Ratio on Core 
Pressure 
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breeding ratio, U 233 critical mass, hydrogen pressure, and 
doubling time. Heat transfer and system analysis studies esti- 
mated the bulk average cavity temperature for the reactor to 
be 3000°K to 4000OK; thermal cross sections and pressures were 
calculated for the case of 4000°K. At this temperature and for 

pressures above about 300 atmospheres dissociation of H ? is 

30 31 

not large; the mole fraction of H 2 is greater than 90% * . 

An exact pressure calculation would be iterative based on the 
H 2 mole fraction, but pressures given here assume an average 
H 2 mole fraction, but pressures given here assume an average 
H 2 mole fraction of 92% at 4000°K. For the doubling time 
calculations a power level of 1000 MW(t) was assumed. A higher 
power level shortens the doubling time proportionally. 

A detailed breakdown of the critical composition for one 
case of the parametric study is shown below: 

H/U = 100/1 , C/Th = 4/1 


Material 

Atom Density 

Mass 

1)^33 

1. 9499x1 0l9cnr 3 

494 kg 

H 

1 . 9499x1 O^ 1 

212 

C 

8.0 xlO 22 

76236 

Th 

2.4xl0 21 

131030 


47 


Some of the conclusions, mostly obvious, within the range 
of this parametric study are noted below: 

1. Breeding ratio decreases as C/Th ratio increases 

2. Breeding ratio appears to be maximum at H/U= 100/1. 

3. Critical mass decreases as C/Th or H/U ratios increase. 

4. Hydrogen pressure increases as C/Th ratio decreases and 

as H/U ratio increases. 

5. Doubling time increases as C/Th ratio increases and 

as H/U ratio decreases. 

The first conclusion can be explained by noting that, as C/Th 
ratio increases, the amount of fertile material decreases 
hence lowering the breeding ratio. The second conclusion is 
essentially an observation but one may note that below a H/U 
ratio of 60/1 the cavity is undermoderated. Above 140/1 the 
effects are much more subtle; increased hydrogen absorption or 
the U resonance may control here, but one desires the lowest 
feasible H/U ratio to yield lower pressure. The decrease in 
critical mass noted in the third conclusion is due to the 
increasing amount of light atoms which soften the spectrum 
toward the large thermal fission cross sections of U^3 # 

Hydrogen pressure is of course expected to increase as hydrogen 
concentration increases, but this also occurs as the C/Th ratio 
decreases. This is because a higher critical mass is required as 
C/Th decreases and, hence, for a given H/U ratio, the hydrogen 
concentration also increases. Variations in the doubling time are 
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due to combined variations in breeding ratio and critical mass, 
with the decreasing mass as H/U increases yielding the strongest 
influence. 

With the trends noted it is difficult to pick an optimal 
configuration. This is because variation in a single parameter 
helps one point but hinders another. For instance, one may 
obtain a lower critical mass by increasing the C/Th ratio but 
this yields a lower breeding ratio also. Or the critical mass 
could be lowered by increasing the H/U ratio, but this, in turn, 
increases the pressure. Of all the cases presented, it can be 
said that critical loadings are within reason; however, pressures 
appear high. 

For the breeder concept one must assess the breeding ratios 
and doubling times. The breeding ratios are low compared to 
that for a fast breeder reactor, but they appear to be reason- 
able for a thermal breeder and yield some reasonable doubling 
times. Breeding ratios near 1.1 with attendant doubling times 
of about 10 years for 1000 MW operation should be quite satisfac- 
tory although engineering details and structure materials will 

probably affect them. In comparison to the molten salt thermal 

32 

breeder with breeding ratios in the 1.05 - 1.07 range, this 
study would show the gaseous core thermal breeder the more 
favorable. For the case of extraterrestrial use, if one merely 
wishes to compensate for process losses, a breeding ratio of 1.1 
should be much more than sufficient. 
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In order to complement the above parametric study efforts 
were turned to examine areas which might improve reactor breeding 
ratio and decrease cavity pressure. One attempt at increasing 
the breeding ratio was by introducing fertile thorium into the 
cavity. For the case of H/U of 100/1 and C/Th of 4/1, thorium 
atoms were added to the cavity in amounts twice, equal and half the 
U 23 3 atom concentration. All three cases resulted in approxi- 
mately a 1 % increase in breeding ratio. Additional cases with 
thorium in the cavity were not examined because of complications 
it would impose on the MHD device. 

The thickness of the blanket region was also examined to 
see if higher breeding ratios could be obtained. Additional 
thicknesses of 20, 40, 75, and 200 cm were examined and the effect 
on breeding ratio is shown in Figure 17. One notes that the 
50 cm thickness used in the study appears to be very near optimum. 
The smallest thickness feasible is desired here to yield lower 
total reactor weights. 

The most obvious method of lowering cavity pressure would 
be to increase cavity size, so a case with a cavity radius of 350 
cm was examined. For a H/U ratio of 100/1 and C/Th ratio of 
4/1, the cavity pressure was reduced from 553 atmospheres to 406 
atmospheres. The breeding ratio also went up 1% due to the 
larger blanket volume. However, the critical mass doubled 
(494 kg to 996 kg). 

One concept which could increase the breeding ratio and also 
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Figure 17. Effect of Blanket Thickness on Breeding Ratio 
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lower cavity pressure is one in which an inner annulus of the 
blanket contains fuel (U 222 ). By placing fissile material 
in the blanket the neutron flux should increase, therefore 
yielding more fertile absorptions, and also reducing the fuel 
required in the cavity for criticality, hence reducing pressure 
for a given H/U ratio. Results for a configuration with fuel 
in the inner 20 cm of a 50 cm blanket region for varying 
quantities of U 222 are shown in Table 5 and Figure 18. One can 
note one disadvantage to this concept, which is that although 
breeding ratio improves with only small amounts of fuel in the 
blanket, pressures are not significantly lower until very large 
amounts of fuel are present in the blanket. 

The use of deuterium as moderator/coolant in place of 

hydrogen could also have the potential of increasing the reactor 

breeding ratio due to decreased absorption (og/og = 1/660 0 

2200 m/s). Since deuterium is not as good a moderator as 

hydrogen higher critical masses would be expected, though. MACH 

233 

estimates for cases with various D/U ratios revealed that U 
masses of 2500 to 5000 kg would be required yielding pressures 
greater than 1500 atmospheres with no case having a breeding 
ratio higher than a comparable hydrogen moderated case. Although 
with deuterium there is essentially no absorption in the moderator, 
the critical mass increases such as to more than offset that loss 
by increased absorption losses in the fuel itself. 



U 233 Atom Ratio 
Blanket/Cavity 

u 233 Cavity 
(kg) 

U 233 Tatal 
(kg) 

Breedi ng 
Rati o 

Fission Ratio 
Blanket/Cavity 

Pressure 

(atm) 

0 

494 

494 

1.100 

.02 

(a) (b) 

553 401 

.5 

489 

553 

1.127 

.04 

552 400 

1 

483 

608 

1.126 

.06 

542 392 

2 

472 

717 

1.125 

.09 

531 384 

4 

468 

919 

1.115 

.15 

529 383 

10 

416 

1496 

1.114 

.26 

465 337 


(H/U-100/1 in cavity, C/Th=4/1; 4000O|<; 250 cm cavity radius; inner blanket 20 cm, C+Th+U; outer blanket 30 
cm, C+Th). 

(a) Hydrogen partial pressure at 4000 °k, H 2 mo le fraction 0.92. 

(b) Hydrogen partial pressure at 3000°K, H 2 mo le fraction 0.99. 


Table 5 Nuclear Data for Gaseous Core Breeder with Fuel in Blanket 


f 



Mass (kg) 



Atom Ratio, 1^233 in blanket / N,j 233 in cavity 

Figure 18. Effect on Critical Mass of Adding U 233 to the 
Blanket Region 
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During the period of these studies a re-evaluation of the operating cond- 
itions of the reactor pointed out that the bulk average cavity temperature 
should be about 3000°K rather than 4000°K. The reasoning for this change is 
that MHD requirements are met with a maximum temperature of about 4000 “K and 
hence the bulk average cavity temperature should be lower. The impact of this 
temperature reduction of the nuclear analysis was shown by a slight decrease in 
critical mass, increased breeding ratio, and, of course, lower pressure. A 
comparison of data for the two temperatures is shown below: 

H/U= 1 00/ 1 , C/TH=4/1 

p-J-3 

Temperature U Mass Breeding Ratio Pressure 

4000° K 494 kg 1.100, 553 atm. 

3000°k 491 kg 1.121 399 atm. 

The primary reason for the lower critical mass and higher breeding ratio is the 
shifting of the thermal neutron spectrum to larger cross section values in the 
"1/v" range. The pressure decrease is essentially linear with temperature, 
but the ^mole fraction at 3000°K also increases to 99%. 

In order to obtain more realistic results additional overall systems 
implications must be integrated into the computations. One important aspect 
is the influence of structural material on the gaseous core breeder reactor. 

In order to estimate such an effect computations were made with molybdenum 
homogeneously mixed in the blanket region. The structural requirements of 
the gaseous core breeder have not been studied, but it is not expected that 
a great deal of structure in neutronically important regions is required. 
However, the following results for cases with structure are shown below with 
the data for no structure: 

H/IM00/1, C/Th=4/1 , 3000 °k 

Atom Percent Mo in Blanket u333^ ass Breeding Ratio Pressure 

0 491 kg 1.121 399 atm. 

0.2 (or 4%, enriched) 493 kg 


1.108 


400 atm. 
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From the above data one sees that for these quantities of structural material 
the breeding ratio is still in the same range as the molten salt breeder 
mentioned previously. One should also note that the absorption loses in Mo 
structure can be reduced by isotopic enrichment in Mo^® and Mo ^0 as notec j j n 
Reference 24. In that case the above results could be equivalent to much 
larger percentages of enriched Mo. 
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GASEOUS CORE NON-BREEDER CALCULATIONS 

This section contains the results of a parametric study of gaseous 
core reactor concepts where breeding of additional fuel is not the primary 
purpose. For extraterrestrial purposes a gaseous core reactor could be 
designed with a sufficiently long life to accommodate many applications. 

The non-breeder reactor may prove to be more favorable than the breeder for 
many applications due to its simplicity and lower total system weight. 

Nuclear calculations have been performed for a range of cavity sizes, 
reflector thickness, and fuels. Cavity radii of 60,80,100, and 150 cm 
have been examined for both 30 cm and 50 cm thick reflectors of beryllium 
oxide (BeO). U^3 f ue -| anc | 35 f ue ] 0 f three different enrichments (.98, 

.93,. 50) have been examined for the various geometries. The bulk average 
cavity temperature is assumed to be 300 °k and pressures are calculated for 
an Hg mole fraction of 99% . Helium would be the more likely coolant for the 
non-breeder but hydrogen was used for expedience. 

Table 6 gives the critical masses for the various cases. The critical 
masses are also depicted in Figure 19 and hydrogen pressures are shown in 
Figure 20. The full matrix of geometric cases was not calculated, rather 
the more likely combinations of cavity radius and reflector thickness were 
examined. Only two cases for the 150 cm radius cavity were examined. 

The results of this parametric study show that the U 2 ^3 f ue -| ec j con _ 

figurations are the most attractive based on critical mass and pressure. 

235 

Critical masses for the U cases are not excessive, but the smaller sizes 
have rather high pressures. As noted above, helium would be the preferred 
working fluid for the gaseous core reactor in conjunction with an MHD device 
and helium would have less absorption than hydrogen. By performing computations 
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TABLE 6 GASEOUS CORE REACTOR CRITICAL MASSES (kg) 
FOR VARIOUS SIZES OF A BeO REFLECTED REACTOR 
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CAVITY RADIUS (cm) CAVITY RADIUS (cm) 

BeO Reflector Thickness = 30 cm 

50 cm 

Average Cavity Temperature = 3000®K 
Hydrogen/Uranium Ratio = 100 



CAVITY RADIUS (cm) CAVITY RADIUS (cm) 

Figure 20. Reactor Pressure for Various Sizes of a BeO Reflected Reactor. 
Average Core Temperature 3000°K. 
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without hydrogen it was found that critical mass and pressure were not overly 
sensitive to the hydrogen as an absorber or moderator. Mass and pressure 
decreased about 5% for the case of no hydrogen. A helium cooled configuration 
should fall between the limits of hydrogen and no hydrogen. 
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Thickness of heat exchanger plate (m) 

MHD duct inlet area (m 2 ) 

Separator inlet area (m 2 ) 

Separator exit area (m 2 ) 

MHD duct i-th segment exit area (m 2 ) 

Frontal area of heat exchanger (m 2 ) 

Free-flow area of heat exchanger (m 2 ) 

* 

Heat transfer area of heat exchanger on hot side (m 2 ) 

Heat transfer area of heat exchanger on cooler side (m 2 ) 

MHD duct i-th segment inlet area (m 2 ) 

Space radiator area (m 2 ) 

Plate spacing of heat exchanger on hot side (m) 

Plate spacing of heat exchanger on cooler side (m) 

Magnetic field strength (tesla) 

Capacity ratio of cool side of heat exchanger }cal/sec-°K) 
Capacity ratio of hot side of heat exchanger (cal/sec-°K) 

Heat capacity at constant presure (cal/gr^K) 

Gas electrical conductivity in i-th segment of MHD duct (mhos/m) 
Distance between two electrods of MHD duct (m) 

Heat transfer effectiveness of heat exchanger 
Friction factor 

Gas flows rate at exit of reactor (kg/sec) 
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Fj Gas flow rate at inlet of MHD duct (kg/sec) 

F 2 Gas flow rate at exit of MHO duct (kg/sec) 

F 3 Gas flow rate at exit of seperator (kg/sec) 

F 4 Gas flow rate at inlet of gas turbine (kg/sec) 

F 5 Gas flow rate at exit of gas turbine (kg/sec) 

F Gas flow rate at exit of first stage compressor which flows into 

8 mixing tank in MODE II (kg/sec) 

Fg Gas flow rate in intercooler (kg/sec) 

F s Liquid sodium flow rate (kg/sec) 

G x Flow-stream mass velocity on hot side of heat exchanger (kg/sec-m 2 ) 

G 2 Flow-stream mass velocity on cool side of heat exchanger (kg/sec-m 2 ) 

h x Convection heat transfer coefficient on hot side of heat exchanger 
(cal/sec-cm 2 -°K) 

n 2 Convection heat transfer coefficient on cool side of heat exchanger 
(cal /sec-cm 2 -°K) 

H Enthalpy of gas (cal/gr) 

H q Enthalpy of gas at exit of reactor (cal/gr) 

Hi Enthalpy of gas at inlet of MHD duct (cal/gr) 

H 2 Enthalpy of gas at exit of MHD duct (cal/gr) 

H 3 Enthalpy of gas at exit of seperator (cal/gr) 

H 4 Enthalpy of gas at inlet of gas turbine (cal/gr) 

H s Enthalpy of gas at exit of gas turbine (cal/gr) 

H g Enthalpy of gas at inlet of first intercooler (cal/gr) 

H ? Enthalpy of gas at exit of first intercooler (cal/gr) 

H 18 Enthalpy of gas at exit of preheater on cool side (cal/gr) 
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H 19 Enthalpy of gas inlet of heat regenerator on hot side (cal/gr) 

H 20 Enthalpy of gas at inlet of reactor 

H m Enthalpy of gas at exit of i-th segment of MHD duct (cal/gr) 

H Stagnation enthalpy (cal/gr) 

k Thermal conductivity (cal/sec-cm-°K) 

K MHD loading factor 
KE Kinetic energy (MW) 

AL.j Length of i-th segment of MHD duct (m) 
m for thin sheet fins 

M x Mach number at inlet of MHD duct 
M 2 Mach number at exit of MHD duct 

M 3 Mach number at exit of seperator 

n Number of segments in MHD duct or number of passes in gas to liquid 
metal heat exchanger 

Np r Prandtl number 

N^ Reynolds number 

P Pressure (atm) 

P Q Reactor cavity pressure (atm) 

Pi Pressure at MHD inlet (atm) 

P 2 Pressure at MHD exit (atm) 

P 3 Pressure at exit of seperator (atm) 

P 4 Pressure at exit of heat regenerator on hot side (atm) 

P 5 Pressure at exit of gas turbine (atm) 

P 6 Pressure at inlet of first intercooler (atm) 
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P7 Pressure at inlet of first stage compressor (atm) 

P 8 Pressure at exit of first stage compressor (atm) 

P 9 Pressure at inlet of second stage compressor (atm) 
P 10 Pressure at exit of second stage compressor (atm) 
P n Pressure at inlet of third stage compressor (atm) 

P 12 Pressure at exit of third stage compressor (atm) 

AP Fractional pressure drop (atm) 

P m Pressure at inlet of i-th segment of MHD duct (atm) 

P m _ Average pressure in i-th segment of MHD duct (atm) 

PR 1 Expansion ratio in MHD duct 

PR 2 Expansion ratio in gas turbine 

PR 3 Compression ratio in each compressor 

PR Expansion ratio in each segment of MHD duct 

m 

Q Reactor power (MW) 

Q. Fraction of reactor power generated in blanket (MW) 
b 

Q c Fraction of reactor power generated in core (MW) 

Q Total heat rejected by each compressor unit (MW) 

K 

Q Heat rejected by first intercooler (MW) 

K i 

Q d Heat rejected by second intercooler (MW) 
k 2 

Q d Heat rejected by third intercooler (MW) 

k 3 

R 3 Radius of heat regenerator (m) 

R 2 Radius of preheater (m) 

Temperature (°K) 


T 
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T q Gas temperature at exit of reactor ( T<) 

T 1 Gas temperature at inlet of MHD duct ( K) 

T 2 Gas temperature at exit of MHD duct {%) 

T 3 Gas temperature at exit of seperator ( ‘K) 

Gas temperature at exit of heat regenerator on hot side ( C K) 

T s Gas temperature at inlet of gas turbine (°K) 

T e Gas temperature at inlet of first intercooler (°K) 

D 

T ? Gas temperature at exit of first intercooler (°K) 

T Gas temperature at inlet of second intercooler (.°K) 

T q Gas temperature at exit of second intercooler (°K) 

T J0 Gas temperature at inlet of third intercooler (°K) 

T Gas temperature at exit of third intercooler ( 0 K) 

T Gas temperature at exit of third compressor (°K) 

T Gas temperature at inlet of mixing tank of cool gas in MODE II (°K) 

T ig Gas temperature at inlet of heat regenerator on cool side (°K) 

T 2q Gas temperature at inlet of reactor (°K) 

Tm. Gas temperature at inlet of i-th segment of MHD duct (°K) 

Average temperature of i-th segment of MHD duct (°K) 

T Temperature of space radiator (°K) 

K 

T r Temperature of liquid sodium at inlet of intercooler (°K) 
i 

T d Temperature of liquid sodium at exit of first intercooler (°K) 

KG ^ 

T d Temperature of liquid sodium at exit of second intercooler (°K) 

T d Temperature of liquid sodium at exit of third intercooler (°K) 

Ke 3 
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T Number of heat transfer units 
u 

T . Heat exchanger hot side inlet temperature (°K) 

1 ,in 

T Heat exchanger hot side exit temperature (°K) 

J ,out 

T , Heat exchanger cool side inlet temperature (°K) 

2 ,in 

T Heat exchanger cool side exit temperature (°K) 

2, out 

AT Temperature difference (°K) 

U Velocity (m/sec) 

U, Overall heat transfer coefficient (Kcal/sec-m 2 -°K) 
h 

V Specific volume (m 3 /kg) 

W Cj Power required by first stage compressor (MW) 

W C2 Power required by second stage compressor (MW) 

W C3 Power required by third stage compressor (MW) 

W Total power output of MHD duct (MW) 

MHD 

W.,,, MHD electrical power output (MW) 

MHD q 

6 

W g Steam turbine power output (MW) 

Wt Gas turbine power output (MW) 

Wax Power demand for auxiliary components (MW) 

X Width of intercooler heat exchanger (M) 

Y Height of intercooler heat exchanger (M) 

Z Length of heat exchanger in direction of gas flow (M) 
a Stefen-Bol tzman constant = 5.67 x 10“ 12 (watts/cm 2 -°K 4 ) 

a Ratio of free-flow to frontal area of hot side of heat exchanger 

1 

a Ratio of free-flow to fontal area of cool side of heat exchanger 

2 
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Ratio of total transfer area of hot side of heat exchanger to 
volume between plates of that side (m 2 /m 3 ) 

$ 2 Ratio of total transfer area of cool side of heat exchanger to 
volume between plates of that side (m 2 /m 3 ) 

Y Ratio of heat capacities 

Y h Hydraulic radius of hot side of heat exchanger (cm) 

1 

Y^ Hydraulic radius of cool side of heat exchanger (cm) 

2 

A Ratio of fin area to total area 
<5 Thickness of fin (cm) 

n Efficiency 

n c Compressor efficiency 

n s Steam turbine efficiency 
n t Gas turbine efficiency 

n sp Cycle thermal efficiency for space power plant 

n G Cycle thermal efficiency for ground based power plant using a 
b t steam bottoming cycle 

Fin effectiveness 

n Q Surface effectiveness 

p Density (kg/m 3 ) 

u Viscosity (kg/sec-m) 
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DESCRIPTION OF THE THREE SYSTEMS 

Three different types of advanced nuclear-MHD power plant systems 
were investigated. The three thermodynamic cycles which were studied . 
are shown in Figures 1-3 and 21-23. The working fluid (hydrogen, 
helium, or argon) is heated in the reactor, passes through a nozzle and 
the MHD generator, and then through two separators. If the coaxial 
flow gas core or colloid core reactor is used, the uranium (and possibly 
• c e s i um - ) particles' would be separated from the working fluid and 
returned to the reactor system. The gas exiting the separator passes 
into four identical heat exchanger and compressor units. One fourth of 
the gas flows through each unit. The three basic thermodynamic cycles 
are referred to as MODE I, MODE II and MODE III. These cycles are 
described as follows: 

1) MODE I (as shown in Figures 1 and 21): The gas from the 

separator passes through a gas to gas regenerative heat exchanger, and 
into a gas turbine which is used to drive the compressor. After exiting 
the turbine, the gas is cooled by a gas to gas heat exchanger (preheater) 
and cooled further by a gas to liquid metal heat exchanger, then 
compressed by a three stage compressor with intercoolers between each 
stage. The gas exits the last stage of the compressor at a pressure 
slightly higher than the reactor pressure, is heated in a preheater, and 
heated further in the reactor blanket and regenerator, before being 
returned to the reactor core. Some of this high pressure gas is 
diverted through a cleanup system to remove gaseous fission products. 



Mj tor He) 


preheater 



MODE X 

•"j 

go 

Figure 21. Terrestrial Regenerative Turbine-Compressor (Mode I) MHD Power Plant. 










MODE II 


Figure 22. Terrestrial Turbine-Compressor (Mode II) MHD Power Plant. 














MODE III 


Figure 23. Terrestrial Regenerative Motor-Compressor (Mode III) MHD Power Plant 
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The sodium coolant removes the heat from the intercooler, and is pumped 
to the space radiator or steam generator. The space radiator is 
employed when the system is used for a space power plant. Sodium-steam 
generators would be used for the ground based power plant. 

2) MODE II (as shown in Figures 2 and 22): After the gas exits 

the separators, it flows into mixing tanks instead of heat exchangers. 

Here it is mixed with cool gas which is taken from the first stage compressor. 
This reduces the temperature of the gas entering the turbine to the 
maximum permissible turbine inlet temperature. Since the mass flow rate 
of gas through the turbine and first stage compressor may be several 
times the mass flow of gas exiting the reactor, the turbine delivers more 
power than in MODE I, and the turbine compressor system is considerably 
larger. Typically the turbine power output is significantly greater 
than the compressor power requirement, so the excess power is used to 
drive an electric generator. The rest of the system is the same as 
MODE I. The advantage of MODE II is that the high temperature regenerator 
is eliminated, so the MHD duct exit temperature can be higher without 
exceeding the permissible heat exchanger inlet temperature. 

3) MODE III (as shown in Figures 3 and 23): This is the simplest 

cycle. The gas turbine and preheater are completely eliminated, and an 
electrical motor is used to drive the compressors. The rest of the 
cycle is the same as in MODE I. The advantage of this cycle is that, 
since the turbine is eliminated, there are no moving parts at high 
temperature. However, this cycle is more sensitive to any degredation 
of the MHD generator performance. If MODE III is used for a terrestrial 
power plant with steam bottoming, the electric motors can be replaced by 
steam turbines. 
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MODE I MAIN PROGRAM 

The computer program MODE I calculates the parameters of the 
advanced nuclear MHD power plant cycle illustrated in figures 1 and 
21. Imput data for hydrogen include: 

1) Enthalpy data for hydrogen at pressures of 1, 3, 10, 30, 

100, 300 and 1000 Atm. for temperatures of 300 to 5000 
degrees Kelvin with 100 degree intervels between the data. 

2) Heat capacity data for hydrogen over the same temperature 
and pressure range, except for 1000 Atm. 

3} Values of the specific heat ratio (y) for hydrogen over 
the same temperature and pressure range 

4) Electrical conductivity of hydrogen seeded with 1 atom 
percent cesium. Helium and argon are assumed to behave 
as ideal gases. 

General Discription 

The data for the enthalpy, heat capacity and heat capacity ratio 
corresponding to a particular pressure and temperature are evaluated 
in sub-pragrams "SBH", "SBCP" and "SBGM" respectively. Given pressure 
P and temperature T, to find the corresponding enthalpy H, the 
following statements are used: 
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IAA = T/100 
AA = IAA x 100 
JJ = IAA - 2 
CALL SBH (P, T, H) 

The first two statements are used to truncate the temperature to an 


ntegral 

mul ti pi e 

of 100°K. 

For exampl 

If 

T = 

531.647 


Then 

IAA = 

T/100 

= 5 


AA = 

IAA x 100 

= 500. 


JJ = 

C \1 

1 

c 
< 
1 — 1 

- 3 


The enthalpy data are stored in an array starting with a temperature 
of 300°K, and data are given for each multiple of 100°K from 300°K to 
5000°K. In the example, the values of AA and JJ allow the subprogram 
"SBH" to select the data of the enthalpy array corresponding to 
temperatures of 500°K and 600°K, and a linear interpolation is performed 
to obtian the value of the enthalpy at 531.647°K. The same approach is 
used to interpolate between enthalpy values given for specific pressures in 
the array* 

Sometimes it is necessary to determine the temperature of the gas 
from known pressure and enthalpy values. This can also be done using 
the same subroutine. For example, given pressure P and enthalpy HY, to 
find the corresponding temperature T, the following statements are used: 

T = estimated value 
IAA = T/100 
AA = IAA x 100 


JJ = IAA - 2 
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1 CALL SBH (P,T,HX) 

DH = HX -HY 

IF [ABS (DH) .LT. lo] GO TO 3 
IF (DH .GT. 10) GO TO 2 
T = T + 1.0 
GO TO 1 

2 T = T - 1.0 
GO TO 1 

3 • CONTINUE 

Keeping the pressure constant, and starting with an estimated 
temperature T, the subprogram "SBH" is used to evaluate enthalpy HX. If 
HX is greater than (less than) HY, the temperature T is decreased 
(increased) by 1°K. This continues until HX = HY within 10 cal/gr. 

Then T is the temperature corresponding to pressure P and enthalpy HY. 

In order to evaluate the parameters of the cycle, 5 initial values 
of temperature, T 4 , T 6 , T 12 , T 18 , and T 20 , are chosen. These initial 
values are used in evaluating heat exchanger characteristics and gas 
properties. After the cycle parameters are calculated, the new values 
of temperature are used and the program continues to Iterate until the 
final solution is reached. 

Discription of the Model 

The heat generated in reactor blanket, Q b , is Q x 0.1, where Q is 
the total reactor thermal power in MWt. The heat generated in the 
reactor cavity, Q c is Q - Q b . 
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The flow rate of gas at the exit of the reactor cavity is 


Q r x 2.389 x 10 2 


H 0 - H 20 


kg/sec 


The enthalpy of the hydrogen is calculated at the inlet and exit of 

the reactor cavity to be H and H , respectively, by subroutine "SBH". 

20 o 

The static temperature and pressure at the exit of the nozzle for 
an isentropic process are 


v • 


i + 2 

2 l 




Where the specific heat ratio y of hydrogen is determined by the 
subprogram "SBGM" at the average temperature and pressure in the nozzle- 
The velocity of the gas at the inlet of the MHD duct is 

U 1 -fi R V' M 1 

The kinetic energy of the gas is 


K.E. = 2 x 2,389 x 10 -1+ cal/gr. 

10 percent of cooler gas with enthalpy H enters through the walls 

20 

of the nozzle for film cooling. An additional 10% is assumed to enter 
the MHD duct walls to provide film or transpiration cooling of the 
MHD duct. The total enthalpy of the mixture is 


(H 0 + 0.1 H 20 ) 
_____ 


cal/gr 


The static enthalpy at the inlet of the MHD duct is 


Hi = H - K.E. cal/gr 
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Knowing the pressure P f = P ( and the enthalpy H, one can evaluate 
the corresponding static temperature using subroutine SBH. The 
density of the gas at the inlet of the MHD duct is 
Pjd.Ol X 10 5 ) 

Pi = Kg/m 3 

RT i 

The mass flow rate at the inlet of the MHD duct is 
Fj = F 0 (l •+ 0.1) Kg/sec 


The MHD inlet area is 


A, = 7; 

1 U Oj 


Dividing the MHD duct into 15 segments and assuming the total expansion 

ratio is PR lf the pressure ratio for each segment PR m is taken to be 
i/n 

PR = PR, 
m * 


Where n is number of segments. 

The pressure at the exit of each segment is 

P = P /PR i = 1,2,3. ..15 
m. „ m. m 

i+l i 

P = P 
m x l 

The pressure drop for each segment is 

AP. = P - P i = 1,2,3. ..15 
i m m 

i+1 i 

The exit temperature of the i-th segment with expansion at constant 

velocity is 

/ 

T 

m 

i+i 

Where K is the loading factor, and y is the specific heat ratio 
corresponding to the average temperature and pressure in each segment. 
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The enthalpy H corresponding to P and T is calculated in subroutine 
m. m. m. 

i i i 

SBH. 

Due to the 10 percent film cooling for the MHD duct, the cooler gas 


with temperature T and static enthalpy H' flows into the MHD duct 

20 20 

and the average enthalpy of the mixture is 


H m.j - 


r / <i A . 0.1, 

f, x (i + i x -rs) 


cal/ gr 


i = 1,2,3. .. 15 

The average temperature T m . corresponding to is found in the 
same way as before. 

The average electrical conductivity of the gas, , in each segment 
is found in the subroutine "SBC" by giving the average temperature and 
pressure 

* _ T + T 

T - m. _ m. k 

m. i+l i 

i 

2 

P + P 

'Xi m i+l m i 


assuming the magnetic flux density for each segment is B. Then the 
length of each segment of the MHD duct is 


aP. x 1.01 x 10 5 

i 

AL i = B2 U c (1 - K) 


i = 1 ,2 ,3. . . 15 


The density corresponding to T m and P m is 


P m _ x 1.01 x IQ 5 


i = 1,2, 3... 15 Kg/m 3 
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The inlet area (Aj^) of each segment (i) is 


I . 

i 


F, + (0.1/15)Fj(1 - 1) . ... , 

_ _ i = 1,2,3. ..15 m z 


U 

1 i 


The exit area (A_) of each segment (i) for constant velocity 

. . 1 
expansion is 


E . 

i 




PR m 


(K^-- 1) 


1 = 1,2, 3. ..15 m 2 


The Mach number at the exit of the MHD duct is 


H = 1 

where T 2 = T is the exit temperature of the MHD duct. 
m 16 

Leaving the MHD duct, the gas enters two separators, if the reactor 
is of the coaxial flow gas core or colloid core type. The Mach number 
of the hydrogen is reduced to M =0.1 before it enters the turbine. 

The temperature and pressure at the exit of the separators are 



(i + *ir m 2 2 ) 

y-T Z 7 

(l + M3 ) 


°K 


P 3 



(1 + ^ Hz 2 ) 1 '" 1 
(1 + 


Atm 


Assuming no heat loss and no frictional losses in the nozzle and 
cyclone separators, the decrease in the enthalpy of the hydrogen passing 
through the MHD duct is equal to the electric power produced. The 
total thermal energy in the -HMB- duct is 

W MHD ■ (CF 0 H 0 ♦ 0.1F 0 H 2q + O.IF^,) - (1. lF^)] x(4.187 x IQ 3 ) MWE 
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Where F Q and f 1 are the gas flow rates at the inlet of the nozzle 
and MHD duct respectively, H 0 is the stagnation enthalpy of the gas at 
the nozzle inlet, H 2Q is stagnation enthalpy of the cooler gas which 
flows through the wall of the nozzle and the MHD duct for film cooling, 

is the stagnation enthalpy of the gas which exits the separators, and 

3 

the numerical value 4.187 x 10 converts Kcal/sec into MW. If two 
separators are used at the MHD exit, then the mass flow rate in each 
separator is 

F, 1.1 

F 2 = — Kg/sec 

The gas passes through a diffuser before entering the separator, 
and the velocity is reduced from M 2 to M 3 . 

The gas velocity at the entrance of the separator is 
U 3 = JtRT 3 M 3 m/sec 

The density of hydrogen at the inlet and exit of each separator is 

p„ = Psd.pl x 10 ) Kg/m 3 
T 3 

The inlet area of the separator is 



Each separator has two exits connected with two turbine-compressor 
units. The exit area is 



where F 3 = — ^ i s the mass flow rate in each turbine compressor unit. 

Before the gas enters the turbine it is cooled by a cylindrical 
counter- flow heat regenerator. The inlet temperature of the hot gas 
is T 3 and the inlet temperature of the cooler gas is T lg which is 
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calculated as follows: 

A value T.. is assumed at the beginning of the program. It will be 
replaced by a calculated value at a later stage of the calculation. The 
corresponding enthalpy H 1S is found from subprogram "SBH". Then the 
temperature T ig corresponding to H J9 can be found. 


H = M + Q b (0. 2398 x 10 3 ) 

" 19 n l 8 

‘2 

Given the heat regenerator radius Rj and length Z l9 the average 
viscosity of hydrogen in both sides and y 2 , the average Prantl number 
N and N , specific heat C p and C D (found from subroutine SBCP), 
mass flow rate F 3 , inlet temperature T 3 and T lg and inlet pressure P 3 
and P 19 , one can calculate the exit temperature T 4 and T 20> percentage 
pressure drops aP 3 , aP 19 , the overall heat transfer coefficient and 
heat exchanger effectiveness using subprogram SBRG . 

Given the turbine exit temperature T 5 , the expansion ratio of the 
gas turbine PR 2 for an adiabatic process is 


Y 

|Y-1 


PR 2 = 


and P u 

P — • — ^ 

5 PR„ 


Atm 


where P 4 = P 3 (1 - aP 3 ) Atm 


aP 3 here is the fractional pressure drop through the regenerator. 
The output of each gas turbine is 


W 


t 


F 4 (H 4 - H s ) 4.187 

rm 


MW 


Where H 4 and H 5 are found by subroutine SBH , 4.187 is the 
conversion factor from Kcal/sec to KW and 1000 changes the units 
from KW to MW. 
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The exit gas from the gas turbine enters a heat exchanger to 
preheat the hydrogen coming from the last stage compressor before it 
enters the blanket of the reactor. The subprogram SBRG is used 
as before. The input data are the size of the heat exchanger (radius R 2 
and length Z ), average viscosity of hydrogen u 1 and (over the 
temperature range of interest), average Prantl number N pr and N pr , 
specific heat C and C D , mass flow rate F , inlet temperature T and 
T , and inlet pressure P 5 and P^. We can calculate the exit temperatures 
T and T._ (this calculated value T 0 is substituted for the previous 
estimated value T,. in the subsequent interation) fractional pressure 
drops aP 5 and aP 12 , the overall heat transfer coefficient and the heat 
exchanger effectiveness. 

The exit hydrogen from the hot side of this gas-gas regenerative 

heat exchanger is further cooled by a gas to liquid sodium heat exchanger 

before it enters the first stage of the compressor. It is a rectangular 

four-pass cross-flow heat exchanger as shown in figure 5. Another 

subprogram SBHE is used for this calculation. The input data are: 

the size of the heat exchanger X , Y , Z , the mass flow rate of the 

hydrogen F , the temperature T , the pressure P where P r = P c (1 - AP ) 

and the mass flow rate F and inlet temperature T of the liquid 

s i R t 

sodium. The output data from this subprogram are the overall heat 

transfer coefficient, heat exchanger effectiveness, exit temperature T y 

of the gas, and T on both sides, the pressure drop aP on the gas 
Kg ^ 8 

side and the pressure head on the liquid metal side. The heat removed 
from the gas side is calculated in the main program. The heat removed 
from the gas side is 

Qr " F 3 ( H 6 " H 7^ 

i 


Kcal/sec 
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The compression ratio for the first stage compressor is 
P, 


P*3 =-p 


B 


where 

and 


m SI 


12 “ 


P 2 P 
6 12 


19 


Atm 

Atm 


1 - &P 1 2 

The exit temperature with isentropic compression 

, T ? 

t 8 

pr 3 ^ 


«K 


The isentropic input power needed for the first stage compressor is 
. _ F a (H' - H ? ) 4.187 


1000 


MW 


For a compressor overall efficiency n c - 0.87, the actual i.nput power 
needed for the compressor is 
W = w' /n 

c i c i c 


The enthalpy of the hydrogen at the exit of the compressor becomes 
H g = 1000 W /(4.187F g ) + H ? 

and the corresponding temperature T can be calculated by subroutine SBH. 

8 

Three sodium to gas intercoolers and compressors are used for each 
gas turbine unit. The calculational precedures are the same as described 
previously. 

The total heat removed from each turbine compressor unit is 
Q = Q + Q + Q p Kcal/sec 

r\ r\ ^ K 2 •'2 

where Q and Q„ are the heat removed from the intercoolers before the 
R 2 r 3 

second and third stage compressors. The total mass flow rate of liquid 
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sodium from each turbine-compressor unit is 


F = F + F + F 
s S S S 

S *2 3 


The mixed temperature of liquid sodium at the exit of the three 
intercoolers is 


F s T Re + F s T Re + F s T Re 

1 1 2 2 3 3 


'r " F 

K e s 


The mixed intercooler exit temperature of the liquid sodium equals 
the inlet temperature of the space radiator when it is used for a space 
power plant divided into 10 regions with the same temperature difference 


AT in each region, 
R 


ATr = 


10 


°K 


T R is an input parameter for the program, 
i 

The average temperature in each region is 


j = 1,2,3. ..10 


t r> = T r - AT R ( j - 0.5) °K 
The total radiator area for each turbine compressor unit is 


A = 
R 


x 4.187 x 10 3 
aE 


X- 

3R T 




When this nuclear-MHD conversion cycle is used as a ground based 
power plant, the heat rejected from each stage of the intercooler can 
be used for steam generation. The steam can be used to generate 
additional electric power or to drive the compressor. The work output 
of the steam cycle is 
w s = Q R X n s 


where Q R is the total heat rejected from each turbine compressor unit 



89 


and n s is the thermal efficiency of the steam power plant. n s is a 
function of the steam temperature. 

The cycle thermal efficiency for the space power plant is 

X ■ [ W MHD + 4 * W t - 4 * + w c 2 + w c 3 >] n 

The cycle thermal efficiency for the ground based power plant is 

n » Tw +4x(W+W)-4x(W + W +W)1/Q 

G. L MHD t s c c c J 

t 12 3 
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THE MODE II MAIN PROGRAM 


All the input data are the same as for MODE I. The heat regenerator 
between the cyclone separator and the gas turbines is replaced by a 
mixing tank. 

The enthalpy at the inlet of the gas turbine is 
,, _ F 3 H 3 + F 1 8 H 1 8 

H ' * " f:v— 

where is gas flow rate at the exit of each cyclone separator, F lfl is 
the gas flow rate of cooler gas from the exit of each first stage compressor. 
This gas is preheated and then enters into the mixing tank to cool the 
hot gas before it enters the gas turbine. The total flow rate of gas 
passing through the gas turbine and first stage compressor is 

kg/sec 

The rest of the gas at the exit of the first stage compressor continues 
through the other intercooling and compression stages of the compressor. 

The gas exits the last stage of the compressor, passes through the 
preheater, and then back to the reactor. 


F 4 = F s + F 1 8 
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THE MODE III MAIN PROGRAM 

All the input data are the same as in MODE I. The gas turbines and 
preheater are eliminated. The inlet gas temperature T, of the first 
intercooler is equal to the exit gas temperature of the heat 
regenerator on the hot side T and the inlet gas temperature on the 
cooler side of the regenerator T ig is equal to the exit temperature of 
the last stage compressor T . In begining the cycle evaluation, only 
3 initial values of temperature T , T and T are chosen. 

r 4 12 20 
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SUBROUTINES FOR GAS PROPERTIES 


Subroutine SBH is used to find the enthalpy of hydrogen corresponding 
to given temperature and pressure conditions. Enthalpy data are read 
into the program and stored in a 48 x 7 array. 


(H) 


H 

H 


H 

— H 

1.1 1.2 

l, 7 

H 

i 

1 

n: 

2,1 2,2 

2 , 7 


} 


1 


I 


cal/gr 


■ i i 

H H H 

I., 48,1 48,2 48,7 

The first subscript represents the temperature and the second 

subscript represents the pressure. The range of temperature is from 

300°K to 5000°K. The enthalpy data are given in each 100°K interval 

at 7 different pressure conditions (1, 3, 10, 30, 100, 300, 1000 Atm). 

For example, the element H represents the enthalpy of hydrogen at 

i ' 

temperature 300°K and pressure 1 Atm., the element H represents 

2 j 3 

the enthalpy at temperature 400°K and pressure 10 Atm., etc. Other 
values of H are found by linear interpolation. For example, to find H 
at T - 325°K, P = 2.5 Atm 


H T = T 1,i + - 

U! (325 - 

100 

300) 

cal/gr 

H 

- H 



H = H + - 

P l ,2 

2,2 1,2 

100 “ 

300) 

cal /g r 
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H - H 

H = H y + P 2 -’ - - (2.5 - 1) cal/gr 

where is the enthalpy at T = 300°K P = 1 Atm 

is the enthalpy at T = 400 °K P = 1 Atm 

is the enthalpy at T = 300°K P = 3 Atm 

H is the enthalpy at F = 400°K P = 3 Atm 

2 ,2 

H y is the enthalpy at T - 325°K P = 1 Atm 

Hp is the enthalpy at T = 325°K P = 3 Atm 

H is the enthalpy at T = 325°K P = 2.5 Atm 

The numerical value "100" in the denominator of the first and 
second equations is the temperature interval between two given data 
at the same pressure and the numerical value "2" in the denominator 
of the third equation is the pressure difference between 3 Atm and 
1 Atm. The numerical value "(325 - 300)" is the temperature difference, 
and the value "(2.5 - 1)" is the pressure difference. 

The subroutines SBCP (for heat capacity), SBC (for electrical 
conductivity), and SBGM (for the heat capacity ratio) use the same 
procedure of linear interpolation as subroutine SBH. The data for 
the heat capacity, heat capacity ratio and enthalpy are taken from 

33 15 

Patch , and the electrical conductivity data are taken from Rosa . 
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HEAT EXCHANGER SUBROUTINES 


Subroutine SBRG 

The subroutine SBRG is used to calculate the cylindrical gas to gas 
counterflow type heat exchanger performance. Input data include: 

1) Size of heat exchanger = Radius R, Length Z 

2) Mass flow rate F 

9 

3) Inlet temperature and pressure on both hot and cool sides 

4) Viscosity of hydrogen gas y and y 2 on both sides 

5) Prantl number of hydrogen gas N ,and N , on both sides 

” r i 

6) Average specific heat of hydrogen gas C and C on both sides 

P 2 

Output data include: 

1) Overall heat transfer coefficient, U^ Kcal/sec-m 2 -°K 

2) Number of heat transfer units T^ 

3) Exchanger effectiveness E 

4) Exit temperature on both sides 

All the surface properties of the heat exchanger and the reference 

data used in the subprogram are taken from reference 34. The subscript 

"1" represents the hot side and "2" represents the cool side. The 

surface characteristics are listed in Table 2. The frontal area of the 

heat exchanger is 

A, = ttR 2 
fr 

The ratio of total heat transfer area of one side to total heat exchanger 
volume are 
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Mi 

“ b i + b 2 + 2a 


ft 2 /ft 3 


b & 

a = — ft 2 /ft 3 

b + b + 2a 
1 2 

Where b and b are the plate spacing of both sides, e and s are ratios 
12 12 

of transfer area to volume between plates, a = 0.012 is used in this 
program .is the thickness of the plate. 


The heat transfer areas for both sides are 


\ ' A fr Z “> 

ft 2 

^h ^fr^ a i 

ft 2 


Where Z is the length of the heat exchanger. 
The free-flow areas on both sides are 


A 

A 


ff. 


ff 


2 


“ Vh/fr 
= Vh A fr 


ft 2 

ft 2 


Exchanger flow-stream mass velocities are 
F 

G = -jr-3- lb/ hr-ft 2 

1 A ff 1 


G 2 = lb/hr-ft 2 

ff 2 


The Reynolds numbers are 

4 \ G i 

N = 1 
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N 


4y G 

V 2 

2 


R. 


2 yj 

where is the hydraulic radius and y is viscosity. The relation between 
Reynolds number and heat transfer characteristics can be expressed in 
two approximate equations derived from experimental correlations given 
by Kays. 34 

h 2/3 (- 0 - 1 " 0 -735 log 10 N ) 

— 1 — n = in 1 

GjC p "Pr! 

2/3 . iQ ( 0*0817 - 0.809 log 10 N R J 


0 2 C„ "Pr 


for the range of Reynolds numbers from 300 to 1000. Both C p and N pr 
are input data, the values of G and N p are calculated from the previous 
equations. So the value of the unit conductance for thermal-convection 
heat transfer, h, on both sides is calculated. 

The correlations 34 for the friction factors in the same range of 
Reynolds numbers as before can be expressed as 


f 1 = 10 


f 2 = 10 


(0.88 - 0.87 log 10 N ) 


(0.9283 - 0.9145 log 1Q N R ) 


The fin effectiveness is calculated from 

tanh (m 1 x jt ) 
n f 1 = m"i x Ii 

tanh (m 2 x a 2 ) 


m 2 x i 2 


'2 ‘ 
ou 

where m =--^- and 


k is the heat conductance of the fin 


6 is the thickness of the fin 


C 


n. 

& 
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The surface effectivenesses are 

v = 1 ■ (1 ■ v 

X ' 1 • x * (1 • n f 2 > 

where the x 1 s are the ratio of fin area to total area for both sides. The 
overall coefficient of heat transfer neglecting the very small wall 
resistance is 


1 = 1 + 1 

u i no h T ( A h /^h J n h 

1 2 "l u 2 

The capacity rate is 

C = F C Btu/hr-F 

h g Pj 


C = F C Btu/hr-F 

c 9 P 2 


The number of heat transfer units is 


Tu = 


A U 
hi i 


T 


The relation between exchanger effectiveness E and number of heat 
transfer units can be expressed as 


n n 
7T 


V E i j/i (Tu - Tu i } 

E = > H 

“ TI (Tu. - Tu,) 

1=1 j/i 1 J 

where the E/s corresponding to the Tu^s are known. 



The exit temperatures of the gas for both sides are calculated by 


T e 


C h^ T i .in ~ ^ltQut^ = C c^ 2 .out ~ T ^,in^ 


C min^i,in " T 2,in^ Snin^i.in “ T 2,in^ 

The equation for the. pressure drop (neglecting the entrance and exit 


loss) is given as follows 
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where v. and v Qut are specific volumes of the gas at inlet and outlet 
and v is the mean specific volume. 


Subroutine SBHE 

This subroutine is used to calculate the performance of the rectangular 
4-pass gas to liquid metal cross-flow type intercooler as shown in figure 5. 
The surface characteristics are listed in table 3. The relation between 

Reynolds number and friction factor on the liquid metal side is given 

6 2 

by f 2 = 0.46 N r “ ’ , the relation between Reynolds number and friction 
factor on the gas side can be expressed by 


(0.23 - 0.559 log N d ) 
f x = 10 10 R i 


for the range of Reynolds number from 300 to 1000. The relation between 

/ hi \ 2/3 

Reynolds number and ( gV also can be expressed by a 


correlation in the same range 'bf Reynolds numbers. 


GiC p 


■Pr! 


2 /3 (-0.38 - 0.534 log 10 N R ) 

= 10 1 


These equations are taken from experimental correlations by Kays. 
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The calculational procedure for the. overall heat transfer coefficient, 

the number of heat transfer units and exchanger effectiveness are as 

used in subprogram SBRG , except the exchanger effectiveness is 

modified by the number of passes 
N E 

E ' = 1 + E {N"~ 1) 

where N = 4 is the number of passes. 

E is the exchanger effectiveness of a single pass. The equation to 
calculate exit temperatures and pressure drops are the same as used in 
subroutine "SBRG". 
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RESULTS 


Figure 24 illustrates the effect of reactor exit temperature and . 
space radiator temperature on the overall thermal efficiency of a 
regenerative turbine-compressor (MODE I) power plant system. The 
upper solid curves are for a terrestrial cycle. The dotted lines 
represent the same cycle but with heat rejection from a space radiator. 

The lower solid curves are the total space radiator area. 

As the radiator temperature is decreased, the efficiency of the 
space plant increases but so does the size and weight of the radiator. 

The final choice of heat rejection temperature will depend on an 
economic analysis of the whole system to determine the optimum 
compromise between radiator size and efficiency. The radiator size 
decreases as the reactor temperature is increased due to the increase 
in plant efficiency with reactor temperature. As the efficiency increases, 
more electric power is produced and less heat is rejected. 

Figure 25 illustrates the effect of MHD pressure ratio on plant 
efficiency. For MODE I the efficiency appears to be insensitive to 
pressure ratio above a pressure ratio of about 4. However, low 
pressure ratios result in high regenerator inlet temperatures which 
adversly affect the reliability of the regenerator. Thus higher 
pressure ratios, of 10 or more, are desired to reduce both the 
regenerator temperature and the turbine inlet temperature. The high 
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AVERAGE TEMPERATURE OF HEAT 
REJECTION TO STEAM GENERATOR 
OR RADIATOR (°K) 


MODE I H 2 + 1 % Cs (or equivalent) 


Reactor Power 
Reactor Pressure 
MHD Loading Factor 
MHD Expansion Ratio 


10,000 MW 
200 Atm 
0.8 
15 


Figure 24. MODE I Power Plant Efficiency and Space Radiator Area 
vs. Reactor Exit Temperature and Average Temperature 
of Heat Rejection from Radiator (Space) on to Steam 
Generator (Terrestrial). 
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Figure 25. 


MODE I Plant Thermal Efficiency vs 
(Hydrogen) 


MHD Pressure Ratio 
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space plant efficiency results from the relatively low (750°K) average 
radiator temperature, A radiator temperature increase of 200°K would 
reduce the efficiency by about 10%. 

Figure 26 illustrates the dependence of MHD power output, compressor 
power requirement, turbine power output, mass flow rate of hydrogen and 
overall plant efficiency on the reactor exit temperature for a specific 
MODE I configuration. The mass flow rate drops by more than a factor 
of two as the reactor exit temperature increases from 2500 to 4000°K. 

This results in a corresponding decrease in the compressor work 
required, and an incr ea s e in turbine power. 

Large MHD pressure ratios result in small turbine pressure 
ratios and a high ratio of MHD power to turbine power. The pressure 
ratios can be chosen to make the turbine power equal the compressor 
power required. The plant efficiency (total power output per thermal 
kilowatt) increases with reactor temperature even though the MHD power 
decreases due to the reduced mass flow rate. Lower flow rates also 
result in smaller compressors and turbines. 

Figure 27 presents the efficiency for MODE II terrestrial and space 
power plants and the space radiator area as a function of the heat 
rejection temperature and reactor temperature. For a given set of 
conditions the MODE II configuration is less efficient than MODE I 
since the regenerator has been replaced by irreversible mixing of the 
hot gas from the MHD duct with cooler gas from the compressor. This 
recirculation increases the gas flow through the turbine and thereby 
increases the size and weight of the turbine and first stage compressor. 
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AVERAGE TEMPERATURE OF HEAT REJECTION 
TO STEAM GENERATOR OR RADIATOR (°K) 


MODE II H 

Reactor Power 
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MHD Loading Factor 
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Figure 27. MODE II Power Plant Efficiency and Space Radiator 
Area vs. Reactor Exit Temperature and Average 
Temperature of Heat Rejection from Radiator (Space) 
on to Steam Generator (Terrestrial). 
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The effect of MHD duct presure ratio on overall plant efficiency 
is illustrated by figure 28. The dependence on efficiency of an MHD- 
turbine cycle is insensitive to MHD pressure ratio for values greater 
than 5 since the MHD efficiency is taken to be almost as high as the 
turbine efficiency. However, larger MHD pressure ratios result in 
reduced recirculation of gas from the first stage compressor back 
through the turbine, so the size and weight of the turbine and first 
stage compressor is reduced. 

The effect of reactor exit temperature on MHD power output, net 
plant power output, compressor power and turbine power output for a 
specific MODE II configuration is illustrated by figure 29. 

Figure 30 shows the effect of reactor exit temperature and average 
radiator temperature on a MODE III plant efficiency and radiator area. 

For this particular plant configuration, the efficiency of the space 
power plant drops rapidly as the radiator temperature is increased. For 
a relatively low radiator temperature (750°K), plant efficiency is 
insensitive to MHD pressure ratio (figure 31). The effect of reactor 
temperature on the MODE III plant power output, MHD power, compressor 
power and flow rate are illustrated by figure 32. As expected, the 
major reason for the decrease in power output for the higher radiator 
temperature is the large increase in compressor power required. Increasing 
the reactor temperature decreases the compressor power requirement 
because of the corresponding decrease in mass flow rate. 

The magnetic field strength is held constant over the length of 
the MHD duct. The duct length is calculated by considering the duct to 
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Figure 28. MODE II Plant Thermal Efficiency vs. MHD Pressure Ratio 
(Hydrogen). 
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AVERAGE TEMPERATURE OF HEAT REJECTION TO 
STEAM GENERATOR OR FROM RADIATOR (°K) 


MODE III H 2 + 1% C (or equivalent) 
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Reactor Pressure 200 Atm 

MHD Loading Factor ■ 0.8 

MHD Expansion Ratio 15 


Figure 30. MODE III Power Plant Efficiency and Space Radiator Area vs. 

Reactor Exit Temperature and Average Temperature of Heat 
Rejection from Radiator (Space) on to Steam Generator 
(Terrestrial ) 
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Figure 32. MHD and Net Power Output, Compressor Power Required and 
H 2 Mass Flow Rate vs. Reactor Exit Temperature for 1200°K 
and 800°K Radiator Temperatures, for a Specific MODE III 
Configuration. 
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be devided into 15 segments, each with a pressure ratio equal to the 
total raised to the 1/15 power. The length of each segment is 
calculated, and the total length is the sum of the 15 segment lengths. 
Figure 33 illustrates the relationship between length to exit diameter 
ratio and magnetic field strength for various reactor exit temperatures. 

The MHD duct exit temperature is shown in figure 34 for both 
helium and hydrogen as a function of pressure ratio. Due to its 
higher value of y> the helium temperature drops faster than hydrogen, 
so smaller MHD pressure ratios are used with helium. Figures 35-37 
illustrate the effect of MHD pressure drop on plant efficiency for Modes 
I, II and III operating with helium. Helium would probably be the gas 
used in smaller non-breeder power plants, whereas hydrogen would most 
likely be used for larger breeder reactors. The pressure ratio for maximum 
efficiency is slightly higher for higher reactor temperatures. Figure 
38 presents approximate relations between magnetic field strength and 
L/D ratio. 
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Figure 33. MHD Duct L/D Ratio vs. Magnetic Field Strength for Hydrogen 
+ 1 % Cesium (or Equivalent) and 200°K Nonequil librium 
Ionization. 
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Figure 34. MHD Duct Exit Temperature vs. Pressure Ratio and 
Reactor Exit Temperature. 
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MODE I Helium 

Reactor Power 10,000 MWt 

Reactor Pressure 50 Atift 

MHD Loading Factor 0.8 

Av. Radiator Temp. (Space Plant) 750°K 


Figure 35. MODE I Plant Thermal Efficiency vs. MHD Pressure Ratio 
(Helium) 
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MHD PRESSURE RATIO 


MODE II Helium 

Reactor Power 10,000 MWt 

Reactor Pressure 50 Atifl 

MHD Loading Factor 0.8 

Av. Radiator Temp. (Space) 750°K 


Figure 36. 


MODE II Thermal Efficiency vs. MHO Pressure Ratio (Helium) 
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MODE III Helium 

Reactor Power 10,000 MWt 

Reactor Pressure 50 Atm 

MHD Loading Factor 0.8 

Av. Radiator Temp. (Space Plant) 750°K 


Figure 37. MODE III Plant Thermal Efficiency vs. MHD Pressure Ratio 
(Helium). 
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MODE III 
Reactor Power 
Reactor Pressure 
MHD Loading Factor 
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Figure 38. MHD Duct L/D Ratio vs. Magnetic Field Strength for 
Helium + 0.451! Cesium (or equivalent) and 200°K 
Nonequillibrium Ionization. 
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the space radiator Preceding page blank 


Assuming an emissivfty of 0.9, the total area required for SNPS 
radiators was calculated. An SNPS producing 13,000 MW of electrical 
power at an overall thermal efficiency of 56% would reject about 
10,000 MW of heat, as shown by the upper curve on figure 39. If the 
radiator, as shown in figures 6 and 7, has a length three times its 
width, then the width of the radiator base may be calculated as a 
function of radiator temperature. There are two such radiators used 
with the power plant. As seen in figure 40, the base of a 750°K 
radiator for a 13,000 MWe SNPS would measure 200 meters, and its length 
would be 600 meters. At 1000°K, these dimensions are cut almost in 
half. 




Emissivity = 0.9 


Heat Rejected 
= 10,000 MW 


1000 MW 


100 MW 
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I MW 


Radiator Temperature ( K) 


Figure 39. Radiator Area Required vs. Average Heat Rejection 
' Temperature. 





Length of Radiator Base (meters) 



Figure 40. Width of Radiator Base vs. Average Temperature 

(Length of Radiator = 3 x Width, 2 Radiators used). 
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POWER TRANSMISSION TO EARTH 


A detailed study of the microwave transmission of power from an 
orbiting power station to earth has recently been reported by the 
Grumman Aerospace Corporation 35 and Raytheon 22 . The Grumman-Raytheon 
study considered a system which would transmit 13,000 MW from synchronous 
orbit to provide 10,000 MW of electrical power from the receiving 
antenna on the ground. The transmitting antenna is proposed to be 1 km 
in diameter, and converts high voltage d.c. electric power into a 
3,000 MHz microwave beam with an efficiency of about 90%. Heat produced 
by dissipative power losses in the antenna is radiated to space by 
cooling fins. Atmospheric attenuation of the beam would vary from less 
than 2% on a clear day to about 7% under worst weather conditions. 

This beam would be intercepted at the ground by a rectifying antenna, 
called a rectenna. Schottky barrier diodes uniformly distributed 
throughout the antenna structure provide rectification so that the 
output from the antenna is high voltage direct current. On the basis 
of experiments performed to date, the projected conversion efficiency 
of the receiving antenna would lie in the range of 85 to 90%. Thus, 
depending on weather conditions and the rectenna efficiency, the overall 
transmission efficiency could vary from 70% to 80%. The diameter 
proposed for the receiving antenna to intercept 90% of the power in the 
beam is 6.8 km (4.3 miles) if the rectenna is located at the equator. 
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At higher latitudes the rectenna could be ellipsoidal with a minor axis 
of 6.9 km and a major axis of 6.9/cose km where 0 is the latitude in 
degrees. For example, at a latitude of 40 degrees, the major axis of 
the elliptical rectenna would be 9 km. 

The power density of the microwave beam arriving at the rectenna 
has a gaussian profile, dropping from a maximum intensity of 81 mW/cm 2 
at the center to 8.1 mW/cm 2 at the edge. At a distance of twice the 
antenna radius the power density is 0.009 mW/cm 2 , and at three times 
the radius the intensity is 8 x 10' 8 mW/m 2 . The radiation protection 
guide for humans, as set in 1966 by the American National Standards 
Institute (USAIC95. 1-1966) is 10 mW/cm 2 for continuous exposure, and 
this standard also applies in western Europe. The limits are higher 
for short term exposure. Studies have shown 35 that occupants of 
aircraft which might accidently fly through the beam would not be harmed. 

An exclusion area surrounding the rectenna could prevent humans 
or animals from receiving any significant exposure. Also, since the 
microwave intensity reaching the ground around and beneath the 
antenna is tolorable, this area could be farmed productively. Although 
the rectenna absorbs 99% of the microwave energy striking it, it stops 
very little sunlight, so the land beneath the antenna can be used for 
the production of food. The heat release due to beam attenation in 
the atmosphere, antenna losses, and microwave heating of the land 
around and beneath the antenna is about 10% of the thermal discharge 
from today's most efficient thermal power plants. 
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WEIGHT ESTIMATES 


Ragsdale 36 8 37 has estimated the weight of a gaseous core reactor 
for rocket propulsion, and for a 22,000 MWt reactor with a 3 meter 
cavity diameter and a 76 cm moderator-reflecter, he arrived at 
a moderator weight of 120,000 lbs. and a pressure ressel weight of 
140,000 lbs., based on a reactor pressure of 1000 Atm., which is an 
upper limit of pressures which might be encountered in gas core 
power reactors. Weight estimates for the hydrogen turbopump range 
from a low of 5000 lbs. to a high of 24,000 lbs. Nuclear calculations 
for the gaseous core breeder reactor, given earlier in this report, 
resulted in a moderator weight of 168,000 pounds and a weight of 
thorium fertile material of 288,000 pounds. The weights of hydrogen 
and fissle uranium are almost negligible by comparison: 466 

pounds for the hydrogen and 1086 pounds for the uranium. If the 
total hydrogen and uranium weight in the plant is four times that 
in the reactor core (two times would probably be more realistic), 
then the total uranium weight would be about 4000 pounds and the 
total hydrogen weight about 2000 pounds. Thus, an upper limit on 
the reactor weight can be arrived at: 
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Moderator 
Pressure Vessel 
Thorium 
Uranium 
Hydrogen 

Other Components 


170.000 lbs 
140,000 , lbs 

288.000 lbs 

4.000 lbs 

2.000 lbs 

75,000 lbs 


(77,000 Kg) 
(64,000 Kg 
(130,000 Kg) 
(-2,000 Kg) 
(-1,000 Kg) 
(34,000 Kg) 


TOTAL WEIGHT 679,000 lbs (308,000 Kg) 


If the SNPS power plant is to produce 13,000 MW of electrical 
power at an efficiency of 50%, the reactor must have a thermal power 
output of 26,000 MWt. The weight of a nuclear reactor is not 
proportional to power output; the percentage increase in weight is 
much less than the percentage increase in power output. However, 
adapting the conservative position that the weight is proportional to 
power output, the projected weight of the 26,000 MWt reactor would 
be 800,000 lbs. (363,000 Kg). 

As a comparison, the 1100 MWt NERVA XE-Prime Engine weighs 

40,000 lbs. (18,000 Kg). Based on this power to weight ratio at 

a power of 1100 MWt, any reasonable extrapolation of NERVA technology 

to 26,000 MWt will yield a reactor weight of less than 800,000 pounds, 

39 

even when allowance is made for breeding. Westinghouse conducted 
an engineering study of the colloid core reactor and arrived at a 
weight of 41,000 lbs. (19,000 Kg) for a 2000 MWt reactor. No attempt 
was made to optimize the weight of the reactor. Thus, this reactor 
with twice the power level would have the same weight as the NERVA. 

A linear projection, using this power to weight ratio, to 26,000 MWt 
would yield a total reactor weight of 533,000 lbs. (242,000 Kg). This 
is a very conservative estimate, even when allowance is made for 
breeding. 
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Thus it is seen that, regardless of the type of reactor used 
(solid core, colloid core or gas core), a conservative estimate of 
the total reactor weight is 800,000 lbs . (363,000 Kg) 

Most of the weight of the MHD generator is the weight of the 
superconducting magnet. Rosa 15 has developed techniques for 
projecting superconducting magnet weights for MHD generators of up 
to 10,000 MWe output. For a field strength of 10 Tesla and a flow 
velocity of 1000 m/sec, the magnet weight would be 11,000 lbs, (5000 Kg) 
for an average electrical conductivity of 100 mho/m, which is typical 
for the SNPS system, or 25,000 lbs. (11000 Kg) for a 20 mho/m 
average conductiviy. Stekly 40 , et al , have projected the specific 
weight of a magnet for a 100 MWe generator to be 106 Kg/MW, which is 
about a factor of four higher than predicted by Rosa's correlations. 

If Rosa's correlations are indeed low by a factor of four, then the 
magnet for a 13,000 MWt SNPS would weigh about 100,000 lbs. Since 
most of the MHD generator weight is associated with the magnet, the 
total weight of the MHD energy conversion system would be less than 
200,000 lbs. Thus, 200,000 lbs, is taken to be a conservative 
estimate of the MHD generator. 

Projections of turbi ne-compressor weights to thousand megawatt 
power levels are difficult to make since large turbines have only 
been built for terrestrial power generation and weight minimization 
was not a major factor. Based on the mass flow rate, temperature, 
pressure and velocity of the hydrogen passing through the turbines, 
the turbine volume is calculated to be about 100m 3 . Similarly, the 
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compressor volume is about 100 m 3 , and the heat exchanger total about 
200 m 3 . Using an average material density within the turbine compressor 
and heat exchanger, of 0.1 gm/cm 3 , the total weight of each turbine- 
compressor-heat exchanger unit would be 88,000 lbs. (40,000 Kg). 

There are four such units with a total weight of 352,000 lbs. (160,000 Kg). 
Thus, including structure and piping, the total weight of the turbine- 
heat-exchanger-compressor system is taken to be 400,000 lbs, (180,000 Kg). 

Ragsdale^ 6 made use of a study by Haller^ 1 to arrive at 
a specific radiator weight of 140 Kg/MW for a large size radiator 
operating at 1100°K. The use of advanced heat-pipe radiators should 
reduce this specific weight considerably, but using the valued of 
140 Kg/MW, the weight of the radiator required to reject 13,000 MW 
of heat is 1,820,000 Kg, or 4,000,000 lbs . 

Other system components include the uranium and thorium 
reprocessing system, the radiactive waste storage and ejection system, 
electric motors (MODE III), various pumps, the control system, and a 
shield to protect delicate electronic components from nuclear 
radiation damage. 

O c 

Based on the weight of the NERVA shield, Ragsdale determined 
that a disk shadow shield for a 22,000 MWt gas core reactor would 
range from 180 to 225 gms/cm 2 . Taking the weight to be 225 gm^m 2 , 
the total weight of a 10 m diameter disk shadow shield for the SNPS 
reactor would be 177,000 Kg (390,000 lbs.). With supporting structure, 
this becomes 400,000 lbs. (182,000 Kg). The total weight of the 
fuel and thorium reprocessing system is difficult to project, and 
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any comparison with today's reprocessing plants is unwarrented, since 
in the case of the colloid and gas core reactor there are no fuel 
elements to fabricate and disassemble. It is believed that such a 
facility for the SNPS would probably have a total weight in the range 
of 1,000,000 to 2,000,000 pounds. The value of 2,000,000 pounds for 
the reprocessing system is used in estimating the total system weight. 
Similarly, a value of one million pounds is assumed for the waste 
deposal system, and the total weight of all other plant components 
including pumps, the control system, and electronics, but not 
including the microwave system, is taken to be 1,200,000 lbs. 

Based on these estimates, the weights of the SNPS components and 
total power plant weight are detailed below: 


SNPS 13,000 MW e WEIGHT ESTIMATES 



x IQ 3 pounds 

x 10 3 Kg 

Nuclear Reactor 

800 

363 

MHD Systems 

200 

91 

Turbi ne-Compressor-Heat 

Exc. 400 

182 

Radiator 

4000 

1820 

Shi eld 

400 

182 

Reprocessing System 

2000 

910 

Waste Disposal 

1000 

450 

Other 

1200 

545 

TOTAL POWER PLANT 

10,000 

4,535 

Microwave Antenna 

9000 

4082 

Additional Thorium* 

1000 

453 

TOTAL SYSTEM WEIGHT 

20,000 

9,070 


♦Additional thorium is provided here to permit up to 40 years of 
reactor operation at 26,000 MWt. 



131 


COST FACTORS 


The major cost difference between the SNPS and similar types of 
terrestrial power plants is the fact that it must be assembled and 
operated in orbit. Grumman 35 did an extensive study of propulsion 
requirements for the SSPS and arrived at a "most likely" cost of 
$100/1 b for transporting the system components to synchronous orbit. 
This cost was based on making use of a reusable space shuttle to 
deliver components to low earth orbit and an ion propulsion system 
for transportation from low earth orbit to synchronous orbit. 

The Grumman study noted that a considerable savings in propulsion 
costs could be effected if the SSPS could be assembled in low earth 
orbit and then, when completed, boosted into synchronous orbit. 

However, it was considered impractical to assemble an SSPS in low 
earth orbit because of various factors relating to the size of the 
solar arrays and the effects of the Van-alien radiation on solar 
cells. These factors would not be important for an SNPS, so the SNPS 
would certainly be assembled in low earth orbit and then be transported 
to synchronous orbit after completion. This should reduce space trans- 
portation costs below the $100/1 b projected for the solar power plant. 
Since the projected total weight for the SNPS power plant system, 
including microwave antenna, is 20 million pounds for a 10,000 MWe 
plant, even if the space transportation cost per pound is taken to be 
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$100, the increase in capital cost of the plant due to space 
transportation is $200/KWe. 

The SSPS study 35 also arrived at a total cost of the microwave 
transmission system of $120 per KWe delivered to the ground, and a total 
cost of the receiving antenna and rectification system of $50 per KWe.. 

These systems would be identical for the SNPS, so the costs should be 
the same. 

The cost of the nuclear power plant is difficult to project, but 
should be comparible to present nuclear plants. Colloid core and 
gaseous core reactors require no fabrication of fuel elements. The 
simplified fuel cycle for these reactors offers the potential of 
considerable cost savings in this area. The high power level (about 
25,000 MWt) of the proposed SNPS improves further the economics of 
on-site fuel reprocessing. On the other hand, the requirement for 
remote operation and maintenance will increase operating costs in 
comparison with similar terrestrial power plants. Thus, with the 
costs of conventional nuclear power plants running $300/KW, the 
projected cost of the SNPS nuclear power plant (exclusive of transportation) 
is taken to be $500/KW, or 5 billion dollars for a single SNPS power 
plant. For a system of this size, however, the final cost per plant 
may be considerably less than 5 billion, especially if a number of 
these plants are built. The total capital cost of the SNPS may be 
broken down as follows: 
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Nuclear-MHD Power Plant $500/KWe 
Microwave Transmission System 120/KWe 
Power Receiving System 50/KWe 
Space Transportation 200/ KWe 


TOTAL SNPS CAPITAL COST $870/KWe 


These data indicate that the capital cost of an SNPS may be well 
under $1000/KWe, assuming that a reusable space shuttle is available 
to place plant components in low earth orbit and an advanced nuclear- 
MHD power plant is developed. These cost projections do not include 
any research or development costs, such as shuttle development. 

An important aspect of SNPS economics is that all societal costs 
are internalized to the power plant system. There are no "hidden" 
costs to society associated with pollution or depletion of non-renewable 
resources. Since the reactor would breed its own fuel from fertile 
thorium, plentiful supplies of fuel would be available for the next 
thousand years or more, which provides plenty of time for the development 
of more exotic energy sources, such as fusion with direct conversion. 



APPLICATIONS 


The nuclear-MHD power plant system which has been described may have 
a number of applications in addition to the SNPS power plant. Such plants 
using compact non-breeder reactors could produce power in the multi- 
megawatt range for a variety of missions. Figure 41 illustrates a 
MODE I plant in use for electric propulsion. Figure 42 shows a MODE II 
plant in space. 

Figure 43 depicts a terrestrial MODE II type power plant using 
a coaxial flow gas core reactor. > Any of the other three types of 
reactors could also be used. Turbine-compressor units are shown on 
opposite sides of the reactor, and a sodium-steam generator is 
depicted behind the unit on the right. Cyclone-type separators 
beneath the reactor core help separate uranium droplets from the 
carrier gas. A MODE III terrestrial power plant is shown in figure 44. 

To the right of the reactor is a motor-driven compressor unit 
attached to a high temperature regenerative heat exchanger. One of 
the magnet coils is shown to the right of the MHD duct. The sodium 
steam generator to the left of the reactor is connected to an adjacent 
hydrogen-sodium heat exchanger. 

The simplified fuel cycle and high efficiency of terrestrial colloid 
core or gas core nuclear-MHD power plants of these types offers the 
potential of a significant reduction in the cost of nuclear power. 






.37 



Figure 43, Terrestrial MODE II Power Plant (steam generator 
shown behind turbine-compressor unit on right, 
fuel separators located under nuclear reactor). 





Figure 44. MODE III Motor- Compress or Terrestrial Power 
Plant (motor-compressor-regenerator on right, 
steam generator shown on left). 
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The safety of these reactors is enchanced by V the continuous removal 
of fission products from the gas and from the recirculating fuel. In 
the event of a major accident, only very small amounts of long-lived 
gaseous fission products could be released. The thermal discharge per 
electrical megawatt from such a plant operating with an 80 % efficiency 
is only 1/4 of the thermal discharge from a 50% efficient plant, and 
only 1/6 the discharge from a 40% efficient plant. Also the efficiency 
of a nuclear MHD-power plant decreases only slightly when the heat 
rejection temperature is raised (such as by switching from wet to dry 
cooling towers) whereas other types of thermal power plants are much 
more strongly affected. Thus, for terrestrial power generation, 
large advanced nuclear-MHD power plants offer the potential of low 
cost power with enchanced safety and greatly reduced environmental 
impact. 
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